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Abstract 

We present a classical probability model appropriate to the description of quantum randomness. 
This tool, that we have called stochastic gauge system, constitutes a contextual scheme in which the 
Kolmogorov probability space depends upon the experimental setup, in accordance with quantum 
mechanics. Therefore, the probability space behaves like a gauge parameter. We discuss the 
technical issues of this theory and apply the concept to classically emulate quantum entangled 
states and even 'super-quantum' systems. We exhibit bipartite examples leading to maximum 
violation of Bell-CHSH inequalities like EPR pairs or exceeding the Tsirelson bound like PR- 
boxes, as well as tripartite cases simulating GHZ or W-states. We address also the question of 
partially correlated systems and multipartite entanglements. In this model, the classical equivalent 
of the entanglement entropy is identified with the Kullback-Leibler divergence. Hence, we propose 
a natural generalisation of this function to multipartite systems, leading to a simple evaluation of 
the degree of entanglement and determining the bounds of maximum entanglement. Finally, we 
obtain a constructive necessary and sufficient condition of multipartite entanglement. 

PACS 03.65.Ud (Entanglement and quantum non-locality) 
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1 Introduction 

While quantum mechanics is one of the pillars of modern physics, its very foundations remain con- 
troversial. A crucial point is that probability in quantum theory obeys different rules than do the 
classical probabilities, as emphasized by Richard Feynman in the 1950's. Later, E. T. Jaynes noted 
that quantum probabilities are striking similar to the Bayesian probabilities. One controversial 
issue concerns 'local realism' highlighted by the well-known EPR paradox Nowadays, it is 

widely accepted that local realism does not hold in microphysics. This 'spooky' property is derived 
from a famous theorem by Bell 0, HJ complemented by Kochen and Specker Q and supported by 
a number of experimental verifications As a result, a huge literature has come to light, aiming 
to reconcile this apparent inconsistency, if not with common sense, at least with a form of logic. 
Currently, a number of different approaches are competing, often at the interface between physics 
and philosophy (see for instance references 043 )■ However, following Jaynes [l4| a nd among 
several authors, [15j-[24j we have warned about a flaw in the proof of Bell's theorem 25J|. In order 
to falsify Bell's claims, we have exhibited a very simple classical counterexample with just three 
dice [26J based on the present model. In our opinion, this proves that the concept of quantum 
non-locality remains actually unfounded. Conversely, if quantum mechanics is really compatible 
with classical local realism, quantum phenomena can be classically emulated, at least to some 
extent. That is the purpose of the present paper. 

We have constructed classical tools, namely, consistent parametric probability distributions, 
that we have called stochastic gauge systems. This model constitutes a contextual probability 
scheme in which the Kolmogorov probability space depends upon the experimental setup, in ac- 
cordance with quantum mechanics. Hence, the recourse to non-locality is not necessary. By 
contrast, the flaw in Bell's theorem is to postulate the existence of an absolute probability space 
whatever the settings, in opposition to both theoretical and experimental evidences. Beyond the 
EPR paradox, we suggest that classical simulation of quantum phenomena may provide some 
insights in several fundamental questions. 

Our underlying picture is the following: We propose that the classical equivalent of a sta- 
tionary quantum system is a spatially extended object in dynamic equilibrium. Therefore, in the 
absence of perturbation, its evolution is completely deterministic and there is no room for ran- 
domness. This analogy suggests that quantum collapse can be compared with a classical break of 
equilibrium. The stability of a classical equilibrium is similar to the resistance to decoherence of a 
quantum state. Based on Bell theorem, it is generally assumed that the notion of entanglement is 
strictly of quantum nature [27, 28]. By contrast, we argue that the classical analogue of quantum 
entanglement is the property of contextual dependence, better underscored when the break of equi- 
librium of a spatially extended object can be triggered from two distant locations separated by a 
boundary. We have found that this boundary behaves like an ignition front where the contextual 
gauge probability trial is performed. Therefore, the contextual entropy should be located on this 
boundary. This resembles the holographic hypothesis suggested by 't Hooft [2^ | and Susskind [30l | 
in black hole physics. 

We will first deal with the technical issues of our gauge probability model before addressing the 
question of classical simulation of quantum entangled systems. We will then describe the classical 
Bell-type model and will exhibit a number of examples with several settings and total correlation 
between two regions. Specially, we will compute exactly some typical bipartite configurations lead- 
ing to maximum violation of Bell or CHSH inequalities. We will next extend the model to general 
locally consistent systems, including non quantum devices like PR-boxes and multipartite systems 
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like GHZ and W tripartite states. We also deal with partial correlation and define the concept 
of classical entanglement entropy. The entropy of a classical entangled system is not extensive 
contrary to the entropy of a separable system. In thermodynamics, this supports the existence 
of long range interactions (3ll. [33|. In this context, the classical entanglement entropy should be 
identified with the Kullback-Leibler divergence. We propose natural generalizations, including the 
definition of classical 'multipartite entanglement entropies' and a concept of 'entanglement scheme' 
to characterize complex multipartite entanglement. 

2 Bipartite Bell type systems 

In classical physics, suppose that a spatially extended object is in equilibrium and that this equi- 
librium can be broken from two distant locations. Equilibrium is a global concept: Break of 
equilibrium from one location signifies break of equilibrium for the whole object and consequently 
break of equilibrium at the second location as well. This means instantaneous correlation at a 
distance between the two locations. Nevertheless, this does not imply instantaneous physical effect 
since the relaxation towards a new equilibrium will require a delay. A trivial picture is the evo- 
lution of a stationary huge soap bubble (possibly in dynamic equilibrium, for instance in uniform 
rotation). We compare this phenomenon with what happens in quantum EPR experiment. For 
this purpose, we identify the spatial extension of the quantum system with the locus of all possible 
regions of measurement. The quantum system in unitary evolution is considered as a stationary 
object in dynamic equilibrium. The measurement of one parameter in one location interrupts the 
unitary evolution or breaks the equilibrium. The final state is a new dynamic equilibrium. We 
claim that the behaviour of similar extended objects in equilibrium may be viewed as a form of 
classical entanglement. 

2.1 Bell-type systems 

Consider a pair of quantum entangled entities {£0, £1} (e.g., particles). The pair, or its global wave 
function, will be compared with a classical extended object stretching out on two space regions TZo 
and 72.x- We will first focus on totally correlated entities (as defined in Sec. 12 . 1 .3f) referred to as 
Bell-type systems. 

In region TZq, we suppose that an observer Oo selects freely a setting uo (e.g., an angle of 
polarization), element of a given set ©o, and measures on £0 a random dichotomic observable xq, 
element of a dyadic set X = {0, 1}, (or a spin so = 2xo — 1) . This measurement breaks the dynamic 
equilibrium of the whole extended object. Similarly, in region 7ti, a second observer 0\ selects 
independently a setting u\ G ©1 and measures on £1 a random dichotomic observable x\ € X, 
(or a spin s\ = 1 — 2xij3 For the sake of simplicity, we will assume that ©o = ©1 = ©• 

The spins are well defined physical observables and thus the interchange of the numbering in 
the pair {£0, £1} interchanges the binary digits '0' and '1'. More generally, we will call bit reversal 
the exchange of '0' and T' in X. Due to symmetry between £0 and £1, physical Bell-type systems 
are often invariant by bit reversal. 

Note that for an ensemble of runs, the binary strings {xo} and {xi} can be viewed as digital 
signals emitted from two distant ports TZo and TZi. 

2.1.1 Contextuality and Bell's inequalities 

As soon as the system {£0, £1} interacts with the settings uo or/and ui, instantaneously, the 
initial equilibrium is broken. Therefore, there is some conditional probability P(a;o; a;i|uo; u\) that 
(Do will observe (after a delay) the outcome xo and Oi the outcome xi. In fact, the function 
P(xo; xi\uo; Mi) (or P(xju) in short) describes the symmetry of the extended object. It will be 
convenient to have a special name for the entries (x|u). 

Definition (Target). We will name global target or simply target each entry (x|u). Let X = 
{(x|u)} be the set of all possible targets. 

For the purpose in hand, we can identify X with a coarse-grained description of the potential 
future equilibrium of the extended object. According to Bell's theorem [3(], Bell's inequalities 
hold, but only when the theorem is valid, i.e., in non contextual systems. By contrast, if the 
conditional probability P(xo; xi\uo; ui) depends upon the settings uo and tti, Bell's inequalities 
can be violated. 

Actually, there is a number of Bell's inequalities. We have shown [25j that basically, these 
inequalities derive from the well known triangle inequalities with respect to a convenient metrics, 

x In Sec. [2] where we deal with only two regions, we use the convention so = Ixq — 1 and s\ = 1 — 1x\ because the 
concept of total correlation is then simpler. We will drop this convention in the other sections. 
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namely the Hamming distanc^ dn({xo}, {xi} ) between the digital signals {xo} and {xi} emitted 
from the two ports TZq and IZi for an ensemble of runs. A similar metrics was proposed by 
Santos 34| using Boolean logic. Define the Hamming divergence for one run as 

1 1 

d(u ,m) = E[dn(x ,xi)] — ^ dH(zo;a:i)P(xo;a:i|no,Mi) ; (1) 

xq=0 xi=0 

where E[.] stands for the expectation value with respect to P(x|u). Clearly, we have 
d(u , ui) = P(l; OK; mi) + P(0; l|u ; Ui) < 1. 

The Hamming divergence is invariant by interchange of the two regions, but not, in general, 
by permutation of the setting^. Therefore, if 81,62 £ ©, in general d(#i,#2) 7^ d(#2,#i). Let 
80,81,62 £ © be three different settings: Then, Bell's inequalities simply read: 

d{6 ,8 2 ) <d{6 ,8 1 ) + d(8 1 ,8 2 ). (2) 

In non contextual systems, we have P(xo; :ei|mo; wi) = Prob(a;o|^o)Prob(:ri|Mi), and the proof can 
be directly derived from an elementary calculation. The original Bell's inequalities make use of 
the expectation value S(uo,wi) = E[so • si]. In non contextual systems we have [25| . 

d(w ,wi) = (1/2)[1 + S(u ,ui)]. (3) 
Thus, starting from the following formulation of the triangle inequality, 

|d(0 o ,0i)-d(0 o ,02)| <d(0i,0 2 ), 
we obtain the well known original Bell's inequality, 

|S(0o,0i)-S(0o,fti)| <1 + S(0i,0a). (4) 

Another popular formulation is the CHSH inequality |3a |. that involves four settings, A,A',B 
and B' , with A,A',B,B' £ ©. In regions IZo (resp. IZi), it is possible to select the settings A or 
A' (resp. B or B'). For a particular situation, the outcomes are respectively a, b, a' and b' with 
a, a', b, b' G {0, 1}. We have, 

< d H (o, 6) + d H (a', 6) + d H (o, 6') - d H (o', b') < 2. 

These inequalities are valid by inspection given that all outcomes can only take the values or 1. 
Therefore, they will hold true for any convex combination of these inequalities. This will be the 
case if it is possible to prepare different experiments with different settings governed by the same 
probability space. As a result, accounting for Eq. |T]), we have 

< d(A, B) + d(A', B) + d(A, B') - d(A' , B') < 2. (5) 

or using Eq. ©, 

\S(A,B) + S(A',B)+S(A,B')-S(A',B')\<2. (6) 

In conclusion, Bell-CHSH inequalities are always valid provided that the target probabilities 
derive from the same underlying probability space, irrespective of the contextual parameters tto and 
Mi. In this classical framework, violation of Bell's inequalities is simply a criterion of contextual 
dependency of the probability space and does not in the least imply instantaneous effect at a 
distance nor violation of local realism. 



2.1.2 Partial measurement, nonsignaling correlations, local consistency and 
causal horizon 

Quantum mechanical collapse may be viewed as a break of equilibrium of the wave function. When 
the particle £1 interacts with the setting m in region IZi, the whole system decoheres and the 
pair of particles is no more entangled. Then, possibly after a delay, the wave function splits into 
two local wave functions. However, in region IZo, this splitting is not observable and the observer 
(Do is not aware of what happens in Hi; he independently selects his own setting uo and proceeds 
to the measurement of xo, while the local marginal probability Prob(a;o|wo) of xo given «o is not 
affected, as long as both the setting u\ and the outcome x\ are ignored. In other words, the local 



2 The Hamming distance |33|| between two binary sequences is the number of places where the two sequences differ. 
Therefore, we have d[j({^o}i {^l}) = 5Z( X © ^1)1 where © denotes addition modulo 2. 

3 For this reason, we use the name 'Hamming divergence' for 'mean Hamming distance d(fco,fci) between regions 7?.q 
and 1Z\ for the settings uo = ko and u\ = k\ respectively'. 



4 



probability Prob(a;o|wo) depends on the knowledge of the observer (Do and is thus an observer- 
dependent concept. This property characterizes quantum partial measurements in two different 
locations, when the global Hilbert space is regarded as a tensor product of two partial subspaces. 
This means that partial measurement is not a way to communicate and is thus compatible with 
space-like separation of the two regions. Adoptin g th e point of view of the observers, we will call 
local consistency this nonsignaling correlation 36-38]. The same property is also encountered in a 
number of situations in relativity, cosmology and black hole physics, where quantum entanglement 
holds between two regions separated by a so called causal horizon, and even in some classical 
laboratory analogues of black holes [39l. l4p||. On the contrary, it is rarely met in usual classical 
physics, e.g., a soap bubble is not locally consistent. Nevertheless, as we shall show in the following 
sections, it is quite easy to design ad hoc classical devices, like dice games [26| . which fulfil this 
criterion. 

Consider a system with K discrete settings, = {9o, 9\, . . . , 9k-i}- When local consistency 
holds, the partial probability in one region, say IZo, does not depends on the setting u\ in the 
second region 1Z\. Let no = 9 k be selected in region IZo we have: 

Viii G © : Prob(x |uo = 9 k ) = P(x ; O|0 fc , m) + P(a:o; l|0 fc ; m) (7) 

and similarly for the second port. When the 2-region system mimics a quantum situation, Eq. (J7J 
describes the partial state in region IZo, obtained by partial trace over the region 72,1. 
Accounting for the normalization relation 

Prob^i = 0|u t = 9 k ) + Prob^i = l\u z = 9 k ) = 1 

for i = 0, 1, Eq.lJTJ provides a number of 2(K — 1) relations between the target probabilities P(x|u) 
for each setting 9k and 2K(K — 1) for all settings (Table [T]). 

For a continuous ensemble of settings, when P(xo; xi\uo; Wi) is continuous and differentiable, 
we have for i = 0, 1, 

dProb(xi\uj) = Q 
dui-i 

2.1.3 Total correlation 

In Bell-type experiment the outcomes xo = xi are identical when the same setting uo = ui = 6 6 
is selected in the two regions. Then: 

P(1;O|0;0) = P(O;1|0;0) = O. (9) 

This provides a number of 2 additional relations between the target probabilities P(xju) for each 
setting 9 and 2K for all settings (Table [TJ. We will call this property total correlation. 

For totally correlated and locally consistent systems, it follows from Eqs. Q, Q that for £ G X 
and 9 G 

Prob(:ro = £\u = 9) = Prob(xi = £|«i = 9) = P(£; £|0; 9). (10) 

Note that 'total correlation' is different from 'maximum entanglement'. We will discuss the 
concept of maximum entanglement in Sec. 13.2.51 In general, totally correlated systems are not 
maximally entangled and we will meet later examples of bipartite maximally entangled systems 
which are not totally correlated (Table [T%b and !23H ). When the system mimics a quantum situa- 
tion, totally correlated systems describe pure states. 

2.2 Classical evolution of a Bell-type system 

As soon as the system interacts with the settings «o and/or ui, the initial equilibrium is broken and 
the system relaxes towards a future equilibrium. Violation of Bell's inequalities proves that a simple 
random trial is unable to describe the contextual correlations P(xo; a;i|uo; u i)- Hence, we address 
the following issue: Is it nevertheless possible to implement the joint conditional probability P(x|u) 
in classical physics? In other words, is it possible to classically simulate a quantum collapse? 

We claim that the answer is definitively yes: We are going to show that any locally consistent 
and totally correlated Bell-type system relaxation can be achieved by cascading (1) propagation 
with finite velocity towards a so-called ignition point, (2) random trial at this ignition point and 
(3) backwards propagation, again with finite velocity. We will next identify the locus of all possible 
ignition points with a causal horizon. 

2.2.1 Ignition states, ignition set and projection function 

In this paragraph, we will construct a convenient outcome set, that we will call ignition set, which 
meets the requirements of contextual systems. We will call its elements ignition states. Let us first 
examine discrete systems. 
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Discrete settings Consider a finite ensemble of K settings, 8q, 0i, .. ., 8k-i- In Sec. 12.1.11 
we have compared the set of targets X = {(x|u)} with a coarse-grained description of the potential 
future equilibrium of the extended object. Now, we will similarly construct a convenient fine- 
grained descriptioi^ as follows: In region IZi define the local target (xi\v,i). When Xi = £ and 
Ui — 8k, we write (xi\ui) = (£|#fc)i- Each local target, e.g. (£,\8k)o in region IZo, can be regarded 
as the union of all (global) targets compatible with xo — £ and uq = 8k ■ 

1 K-l 

(x \uo) = (£\0k)o = [J |J (x ;xi\ua;ui). 

xx=0 ui=0 

Conversely 

(a;o;a;i|uo;Mi) = (%o\uo) n (xi\ui). 
According to Eq. (|10|) . we have Prob(£|#fc) = P(£; (,\8k; Ok)- When total correlation holds, the 
marginal probabilities are equal in the two regions and the local targets (f|0fc)o and are 
identical. Let Ti oc = {(£|#fc)} De the set of local targets (relevant for both regions). The cardinality 
of Tioc is 2K. 

Let 3 = {j} be a set of integers. We aim to construct a sequence of card(3) fine grains. Define 
an application of x 3 — > X = {0, 1} or rather an application LT(fc, j): JO, K — 1] x 3 — > X = {0, 1} 

e = n(fe,j). (ii) 

We will call projection function the function H(k,j). Now, we construct the sequence of card(3) 
fine grains Xj by the following intersection of local targets: 

K-l 

j£3hV"" n ( n c*>i)i*fc)- ( i2 ) 

fe=0 

Therefore, the maximum number of distinct fine grains is 2 K and then 3 = [0, 2 K - 1]. Let 

K-l 

fe=0 

be the binary expansion of an integer j, (0 < j < 2 K — 1) and jk the coefficient of 2 ft . A convenient 
definition of the projection function is 

n(fc,j) ( = f) j fc . (is) 

For reasons explained in Sec. 12.2.31 we will call ignition state each fine grain Xj. Clearly, the 
ignition states included in a given global target (xo; xi|wo; u i) are also included in the local targets 
(xo\uo) and (a?i|ui) and conversely. Therefore, we have 

Xj e (x ; xi\u ; m) Xj G (x \u ) n (xi\m). 

It is also convenient to define the index subset, 

{j\n(k,j)=o- (w) 

Clearly, we have 1(0, k) + 1(1, k) = 3. For simplicity, we will often identify k with 8k, j with Xj 
and use the same symbol, 

ifcA) (d = ; n(e fc ,A,) ( = f) n(fc,j). 

In addition define the index subset, 

J(x,u) < =' ) {j / Xj € (x,u) } = I(a;o,wo) C\l{x\,u{). (15) 
Finally, let A be the set of all ignition states, 

A < = > {Ao, Ai, . . . , A card (3)_i}. 

We will call this ensemble ignition set. While the maximum value of card(3) is 2 K , we shall 
see that only 2A' ignition states are necessary and thus, this set is strongly redundant for large 
K. A more convenient set 3 of cardinality just 2K is the set of 'double-plateau' integers Da- 
defined as the integers of [0, 2 K - 1] whose binary expansion is a chain of identical bits except 
for a maximum of one jump. For instance, D 3 = {OlT, ITT, TTO , TOO , 000, 00T} = {3,7,6,4,0,1} 

4 This designation describes a method but does not fit very well, because the number of useful fine grains is often less 
than the number of coarse grains! 
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and D 4 = {0011,0111,1111,1110, 1100,1000,0000,0001} = {3,7,15,14,12, 8,0,1}. For reasons 
explained in Sec. 12.3.51 this particular order in Hk, beginning with K/2 zeros for K even or 
(K — l)/2 zeros for K odd, will be referred to as the trigonometric order. For r £ {0,2K — 1], this 
order defines a double-plateau function Dk(t), 

D k : [0, 2K - 1J -> JO, 2 K - 11 : r ^ j = D K {r) (16) 

For example, we have -D,j(0) = 3 or 7)4(6) = 0. 

Note that the ignition states Xj have no definite probability for K > 1. Therefore, in general, 
they cannot be regarded as hidden variables. Nevertheless, we will define in the next section (Sec. 
I2.2.2[) a number of gauge probability distributions gfc(Aj) on A, each compatible with P(x, u). For 
instance, for K = 3, only 6 ignition states are necessary: The ignition states can be pictured by 
dice sides, while each probability distribution corresponds to a particular biased die [2(|. 



Number of settings 


K 


Number of targets (x|u) 


AK' Z 


Number of local targets card(Ti oc ) 


2K 


Constraints of normalization 


K* 


Constraints of local consistency 


2K(K-1) 


Constraints of total correlation 


2K 


Degrees of freedom 


K x 


Number of gauge distributions gfc(Aj) 


K 


Rank of the gauge linear system 


2K 



Table 1: Degrees of freedom in general bipartite Bell-type systems with K possible settings. Due to total correlation, 
the K 2 degrees of freedom are shared between the two regions 



Continuous settings The discrete projection function Eq. (|11|) is not appropriate for continu- 
ous settings. We have thus to guess a convenient projection function 11(0, A) : x R — > X = {0, 1}. 
The fine-grained description for continuous settings is next defined by: 

AeRi->A= P| (n(0,A)|0). (17) 
see 

We will use such a projection function in Sec. 12.3.51 More generally, the structure of the ignition 
set should reflect the system symmetry. 

2.2.2 Gauge probability distributions 

For ease of exposition and simplicity, we will focus on finite systems with K £ N different settings. 

We have emphasized that the ignition states Xj have no definite probability. Now, we are going 
to define a number of K parametric probability distributions gk(Xj) on the ignition set A, labelled 
k from to K ~ 1. Of course, each distribution has to be compatible with the target probabilities 
P(x|u). We regard each target as an union of ignition states. Therefore, the probability of each 
target P(x|u) will be identified with the sum of the probabilities of its ignition states. There are 
a number of consistent possibilities and two families of solutions. 

A first family is obtained from region IZo'- Let no = 6 k and let pkj be the unknown probability 
of Xj in the probability distribution labelled k. For each setting u\ G and each outcome pair 
xo,xi, we write one compatibility condition: 

Vui, Vx ,Va;i : 2Z ptj = P(x ; xi\u — 9 k ; ui) (18) 

iSJ(x.u) 

where J(x,u) is defined by Eq. (|15p . From Eq. (|14p . we derive 

J(xo,xi,Uo,ui) = I(xo,uo) n 1(2:1, 

We have clearly 

Vui 6 © :I(a;o,uo) = I(x ,u ) n [ 1(0, ui) U 1(1, mj] , 
and thus, local consistency, Eq.©, is automatically encoded. Similarly, we have 

j(o,iAA) = i(o,0fc) n 1(1, fc ) =0, 
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and thus, total correlation is also secured. 

Accounting for local consistency and total correlation, each parametric probability distribution 
labelled k is defined by a linear system of AK — 2(K — 1) — 2 = 2K independent equations and 
card(3) unknowns. There is generally a set of positive solutions. A second family of solutions could 
be obtained similarly from region TZi, by interchanging the indices and 1. Select one solution 
Pkj and define 

gfe(Aj) ( =" p k j (19) 

Clearly, we have for any Ui: 

i i 

^2sk(Xj) =J2J2 P(zo;a;i|0 fc ;ui) = 1. 

j£3 Oil- 

We will call the function gk(^j) gauge probability distribution. 
2.2.3 Classical collapse 

We aim to describe the random system relaxation by use of the ignition set A = {A} and the K 
probability distributions gfe(A). The process is the following: In each region, the observer selects 
his own setting uo or u\. These settings are transmitted with finite velocity towards a so-called 
ignition point 3. Only one setting uq or ui, labelled 8 lsn i is kept, e.g., the first received setting. 
Now, at 3, we perform a trial in the ignition set A using the only probability distribution gigni(Aj) 
to draw a single ignition state Xj i ni . Let j UQ — II(Mo,Aj, nl ) and j ui — n(Mi,Aji ni ). These 
coefficients, j U(l and j ui , are transmitted with finite velocity towards the end regions TZo and TZi 
respectively, where the final outcomes are xo = j uo — II(uo, Aj igni ) and xi = j ui — H(ui, Xj i ,). 

Conversely, to compute the probability P(x|u) we have to collect all ignition states Xj within 
the target (x|u): 

P(a; ;a;i|uo;ui) = ^ gigw(Aj) (20) 

Aj6(x|u) 

where the label of gi gn j is such that 8i gn i is any one of the two settings, uo or ui, and Xj £ 
(xo;xi\uo;u\) means that xo = Tl(uo,Xj) and x\ — II(mi, Xj). 

Now the proof is a straightforward verification of the identity of Eq. (|18[) and Eq. (|20p . By 
construction, Eq.© is automatically satisfied and the choice of uo or u± for 8i gn i to select the 
distribution gi gn i at the ignition point 3 does not affect the overall probability P(xo;xi\uo;ui). 

Ignition front and causal horizon Note first that the process implies the existence of an 
ignition point located between the two end regions. Then, the locus of all possible ignition points 
draws a form of frontier between TZo and TZi, whereas no transfrontier communication is possible 
between the two regions. Such an ignition front is not recognized in quantum mechanical collapse 
since this process is generally assumed to be superluminal and non local [41| . but it could be light 
cones, light-sheets [42j| or null surfaces. In black holes physics, it looks like the event horizon. We 
will thus identify this ignition front with a causal horizon [4^. 

Gauge probability distributions There are two possible probability distributions for the 
trial at 3, defined by the settings of the two end regions. The real process, if any, is not ob- 
servable and the choice between these two distributions is indifferent. Therefore, the probability 
distributions gfe(A) behave like global gauge parameters: The choice of one particular distribution 
correspond to a gauge selection at 3 and this choice does not affect the observable result. This 
gauge invariance is fundamental. In our opinion, it is the core of the EPR paradox. We will call 
stochastic gauge system the ensemble of the ignition set A = {A} and a number of parametric 
distributions gfc(A). 

Analogy with the holographic hypothesis The parametric trial is performed on the 
ignition front, i.e., on the causal horizon and not in regions TZo and TZi. This meets a suggestion 
by 't Hooft and Susskind in black hole physics, namely, the holographic hypothesis |29, 30, 42], that 
space of quantum states in a region must be associated with the two-dimensional boundary rather 
than the volume of the system. Note that this assumption might be checked by experiments on 
entangled entities more easily than in black holes. 
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Consistency conditions Due to gauge invariance, the parametric probability distributions 
gfc(Aj) are not independent (25| . Suppose that the observers Oo and Ox select respectively the 
settings uo and Ui G 0. Let g UQ and g ui be the gauge distributions associated with uo and ui 
respectively. For each run, the ignition state can be obtained by a parametric trial on A with 
respect to either g„ or g Ul . Then, Eq.(0 implies that for i — 0, 1 and x — or 1, 

Prob(ar 4 |«i) = £ = £ g ni (A^) (21) 

Clearly, these conditions mean that Prob(xi\ui) is well defined, i.e., that gauges distributions 
are compatible with local consistency. Conversely, when the consistency conditions, Eg, (12111 . 
are fulfilled, using Eq. (|20p . the probability distributions gfc(Aj) define a single random system 
P(xo; xi\uo;u\). An alternative formulation with expectation values reads, 

E(xi|«i)= J2 &»o(Aj)= E ( 22 ) 
l,s(i|«o) A^efilui) 

We will give some examples, but beforehand, let us summarize our results so far: 
Summary (Emulation of Bell-type systems). Any Bell-type system can be emulated by a stochastic 
gauge system composed of an auxiliary outcome set A = {Xj}, called ignition set, and a set of 
gauge distributions gt(Xj). One distribution gt is associated with each setting 8k- A setting m 
is freely selected in each end region IZi and transmitted with finite velocity towards a so-called 
ignition point I, where each setting is associated with its gauge distribution. Only one of the two 
distributions, gi gn i, is used atT. This choice is called gauge selection. A single trial is performed 
at I to draw a unique ignition state Xj. This result is transmitted with finite velocity backwards the 
end regions. The final outcomes Xi are next computed by a projection function H(v,i,Xj). Gauge 
selection at T does not affect the probability of the final outcomes. The ignition point is located at 
the boundary of the two regions. Therefore, if the two regions are separated by a causal horizon, 
the gauge probability system is located on the horizon. 

In classical systems, we will call this random process classical collapse. 
2.2.4 General 2-setting totally correlated systems 

Consider first a Bell-type system with just K = 2 settings, 8o and 9i. When accounting for 
local consistency and total correlation, it is easily shown that the general system depends on 
4 parameters, qi,i = 1,2,3,4 (Tabled, subjected to the following constraints: < < 1; 
93 > 9i) 92; 94 > 91, 92; and qi + qi > qz, 94. The linear system Eq, ()18[l admits 2 A' = 4 unknowns 
for 4 independent equations. The gauge distributions gfc(Aj) are given in Table [3] Suppose that 
the observers choose two different settings, e.g., uo = 80 and Ui = 9i- Then, at the ignition point 
3, any distribution of the working set can be selected. The trial with respect to this particular 
distribution gives an ignition state Xj. Finally, since uq — 80, we get xo = or 1 depending on 
whether j is even or odd. Similarly, xi — or 1 according to j < 2 or j > 2. Gauge invariance 
means that all statistical properties of xq and x\ do not depend on the gauge selection, because 
they can be computed directly from "P{xo\x\\uo;ui). The entanglement is detected by the CHSH 
inequality Eq.© onljjf] when qs 7^ q±. Suppose that 54 > 93 and select A = do, A' =81, B = 9i, 
and B' — 9o. Then we have, 

S(0 O , Si) + S(0 O , Oo) + S(0 u 0i) - S(0i, 0o) = -2 - 4(g 4 - 93) < -2. 

When 53 = qi, the two distributions go (j) and gi (j) are identical. The concept of gauge distribution 
degenerates and thus each ignition state Xj has a definite probability, g(j), whatever the settings. 
Therefore, the ignition states can be regarded in this case as random hidden variables located at 
the ignition point. 

2.3 Simulation of the EPR spin experiment 

We are going to classically simulate conventional EPR-B experiments, i.e., Bell-type setups for a 
pair of totally correlated spins with particular probabilities computed from quantum mechanics. 
In 2D-configuration, let be the unit circle U(l) = [0, 2ir]. From quantum mechanics we know 
that for any a , 0^ £ the target probabilities are given by the following equations and Table Q. 

P(0; l\0 a ; 9 b ) = P(l; O\0 a ;8t) = (1/4) [1 - cos(0 a - b )]. (23) 
P(0; 0\9 a ; 8b) = P(l; 1|0„; 9 h ) = (1/4) [1 + cos(<9 - 8 b )]. (24) 

Clearly, the system is locally consistent and totally correlated. In addition the two entities are 

symmetrical and therefore the two families of solutions are identical. 

5 In Sec. 13.2.51 we will define a more efficient witness derived from the entanglement entropy. 



9 



(Xo, X\) 


(6*0, 6*0) 


(<?i A) 


(#o, #i ) 

V u , J- / 


(01, 0o) 


(0 0) 


1 — 01 




Ho 


1 — r/1 


(0,1) 








93 — 9i 


94 — 92 


(1,0) 








93 - 92 


94 - 91 


(1,1) 


9i 


92 


9i + 92 - 93 


9i + 92 - 94 


d(uo,«i) 








293 - 91 - 92 


294 - 91 - 92 


S(u ,wi) 


-1 


-1 


493 - 29! - 2g 2 - 1 


494 - 291 - 292 - 1 



Table 2: General locally consistent totally correlated Bell- type system with two possible settings 0q and 9i depending 
on 4 parameters 91,92,93 and 94. The table gives the conditional probabilities P(xo; x\ |«o; ^1) as well as the Hamming 
divergence d(uo,«i) and the mean product spin S(«o,wi) between binary strings. 



3 


go(j) 


gi(j') 





1 - 93 


1 - 94 


1 


93 - 92 


94 - 92 


2 


93 - 9i 


94 - 9i 


3 


9i + 92 - 93 


9i + 92 - 94 



Table 3: Gauge probability distributions corresponding to Table [2] 

2.3.1 2-setting EPR-B experiment 

Consider first a 2-setting EPR-B experiment: Let 80 and 81 be the polarization angles and 801 = 
9i — 60. This is obtained with gi = <j2 = 1/2 and — 54 = (1/4) (3 — cos #01) in Table [2] Then 
Table [3] reduces to Table [5] The two gauge distributions are identical and the ignition states can 
be viewed as hidden variables with a definite probability. Such 2-setting EPR-pairs are referred 
to as 'Bell's states' when 9oi ~ tt/2. 

2.3.2 3-setting EPR-B experiment 

We now consider an EPR-B system with K — 3 different settings, leading to a violation of Bell's 
inequalities Eq.©. 

The linear system Eq. (|18p admits 2 3 = 8 unknowns and 2K = 6 independent equations: Each 
probability distribution depends on 2 arbitrary parameters and therefore the general solution 
depends on 6 arbitrary parameters. We can take advantage of this degeneracy to cancel the gauge 
probability of two ignition states. It is possible to use the 6 double-plateau indices D3, 

Da = {j} = {3, 7, 6,4,0,1}. 

The linear system admits then only one solution provided that all gauge probabilities ptj are non 
negative. 

Let 80, 9i and 82 be the three settings with < 9o < #1 < 92 < tt. Let 8 a i, — 9b — 8a- The 
solutions of Eq. (|18|) define a working set of effective gauge ignition states given by Table [6] The 
maximum violation of Bell's inequality, Eq.([5]) or @, is obtained, e.g., for 80 = 0, 81 = tt/5 and 
82 = 7r/2. This could be checked by a computer simulation using Table[71 



(uo,ui) -> 


(0a, 0a) 


(0a, 0b) 


{x ,xi) 


(0b, 0b) 


(0b, 0a) 


I 






(0,0) 


1/2 


(l/4)(l + cos0) 


(0,1) 





(1/4) (1 - cos 6>) 


(1,0) 





(1/4) (1 - cos 6>) 


(1,1) 


1/2 


(l/4)(l + cos6») 


d(u ,ui) 





(1/2)(1 - cos6>) 



Table 4: Conditional probability P(xo; xi \uo; ui) of a typical EPR-B system for two different settings 8 a and 65 with 

e = e a - e b . 
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j 


4 X Sn(j) = 4 X gi (?) 

OU W / O-L \J / 





1 + cos 6*oi 


1 


1 — COS #01 


2 


1 — COS #01 


3 


1 + cos 6*oi 



Table 5: Gauge probability distributions for a symmetric EPRB system with 2 settings, Oq and 6\. (We have $oi = 
$1 — $(>)• The two gauge distributions are identical 



3 


4 x go(j) 


4 x gl (j) 


4 x g 2 (j) 





1 + COS #02 


COS #01 + COS #21 


1 + COS #02 


1 


1 — COS #01 


1 — COS #01 


COS #12 — COS #02 


3 


COS #01 — COS #02 


1 — COS #12 


1 — COS #12 


4 


COS #01 — COS #02 


1 — COS #12 


1 — COS #12 


6 


1 — COS #01 


1 — COS #01 


COS #12 — COS #02 


7 


1 + COS #02 


COS #01 + COS #21 


1 + COS #02 



Table 6: Working set of gauge probability distributions using only 6 ignition states (out of 8), j £ D3 = {0, 1, 3, 4, 6, 7} 
for 3 settings. This working set is only valid when all entries are positive. 

2.3.3 4-setting EPR-B experiment 

We now consider an EPR-B system with K — 4 different settings, leading to a violation of CHSH 
inequalities Eq.([S]). We proceed similarly as for 3 settings. 

Let 9o, 61,82 and 83 be the four settings. A violation of CHSH inequality is obtained, e.g., 
for 9o — 0,9i = n/8,62 = tt/4 and 8z = 37r/8. The maximum violation is obtained, e.g., for 
0o = 0, 0i = tt/4, 02 = tt/2 and 83 = 3vr/4. 

The linear system Eq. (|18p admits 2 4 = 16 unknowns and 2K — 8 independent equations: Each 
probability distribution depends on 8 arbitrary parameters and the general solution depends on 
32 arbitrary parameters. We will take advantage of this degeneracy to force to zero the gauge 
probability of 8 ignition states. It is possible to keep the double-plateau indices, 

D 4 = {j} = {3, 7,15,14, 12,8,0,1}. 

Let again 9 a b — 9b — a - The linear system admits then only one solution given in Table [8] which 
turns out to be positive in the range of interest (Table [9]). 

2.3.4 Discrete regular X-setting EPR-B experiment 

Finally consider a discrete EPR-B system with K regular settings 

fcelo,Jfr-i] : k^9 k = ^. 

We proceed similarly as above. The rank of the system is 2K and the relevant ignition states Xj 
are defined by the double-plateau function j — Dk (r) , Eq. (|16[1 , with r £ [0, 2K — 1] . If we use 



3 


go(j') 


gi(j') 


Ss(J) 





0.250 


0.349 


0.250 


1 


0.048 


0.048 


0.147 


3 


0.202 


0.103 


0.103 


4 


0.202 


0.103 


0.103 


6 


0.048 


0.048 


0.147 


7 


0.250 


0.349 


0.250 



Table 7: Gauge probability distribution with 6 ignition states derived from Table [6] for the 3 settings: 9q = 0; 0i = 
7r/5; 82 = tt/2 leading to the maximum violation of Bell's inequality, Eq. ((2jl or QJ. 
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j 


4 x g (j) 


4 x gAj) 


4 x g 2 (?) 


4 x g 3 (j) 


n 


1 "T CUb (73Q 


pAc fl. „ I rric! f)i n 


CUb C720 cuts C723 


-L CUb C73Q 


1 

1 


1 CUb U±Q 


T PAC hi „ 

1 cub t/xo 


PAC fir-.-, PAQ hi ,^ r~. 

CUb (721 CUb (720 


PAC t"r 1 — . 1 PAQ hi r~. 

CUb "31 — CUb "30 


Q 

o 


CUb 171Q CUb t/20 


1 CUb 17 12 


1 CUb c 1 2 


PAC Hnr. PAQ hi .-. -1 

CUb "32 CUb C/31 


7 


COb [720 — cub P30 


rr^c hi - f ' r~i A/„ - 
LUo "21 tUh C/31 


1 PAC hir-.r-. 

1 COb P23 


1 pAO fir-.n 

_1_ CUb P23 


8 


COb #20 — COb 6> 30 


COS 6*21 — COS 6*31 


1 - COb #23 


1 - COS #23 


12 


cob 6>io — cob 6*20 


1 — COS 6*12 


1 — COb #12 


COS #32 — COS #31 


14 


1 — COS #10 


1 — cos # i0 


COS #21 — COS #20 


COS #31 — COS #30 


15 


1 + COS 3O 


COS #10 + COS #13 


COS #20 + COS #23 


1 + COS #30 



Table 8: Working set of gauge probability distributions using only 8 ignition states (out of 16), j £ D4 = 
{0, 1, 3, 7, 8, 12, 14, 15}, for 4 settings. This working set is only valid when all entries are positive. 



j 


go(i) 


gi(j') 


g2(i) 


g3(j) 





0.073 


0.177 


0.177 


0.073 


1 


0.073 


0.073 


0.177 


0.177 


3 


0.177 


0.073 


0.073 


0.177 


7 


0.177 


0.177 


0.073 


0.073 


8 


0.177 


0.177 


0.073 


0.073 


12 


0.177 


0.073 


0.073 


0.177 


14 


0.073 


0.073 


0.177 


0.177 


15 


0.073 


0.177 


0.177 


0.073 



Table 9: Gauge probability distributions derived from Table[8] for the four settings: 80 = 0; 6>i = 7r/4;#2 = tt/2;(?3 = 
37r/4. This setup leads to the maximum violation of the CHSH inequalities, Eq. J6t . 



the same projection function U(k,j), Eq. (|62ll . we have, 

n[fc, D K (r)} = n[0, D K (r - k)] ; g k [D K (r)] = g [D K (r - k)] 

where r — k > is computed modulo 2K. 

An alternative formulation is to define a number of 2K ignition states, 

a new projection function 

n D (fc,r) ( = f) n[fc, D K (r)} 

and then 



The functions, 



5k 0) = gk[D K {r) 



5(r) ( ^g?(r), and H(r) ( => II D (0, r) 



with r £ JO, 2K — 1] are sufficient to describe all distributions since 

gk(r) = g(r - k) and Tl D (k, r) = U(r - k) (25) 
where i — k > is computed modulo 2K. It is a simple exercise to derive the following expression, 

1. f|cos(2r + l)a I KsO (mod 2) 

g(r ) = - sin a X < (26) 
5 W 2 \ I cos 2ra I ^ = 1 (mod 2) V ' 

where a = n/(2K), and, 

n _ Jo if 2r < K (computed modulo 2K) 
1 1 \i 2r > K (computed modulo 2K) 

(where the residues modulo 2K are taken into [0, 2K — 1]). For instance, Table [lOl gives the result 
for K = 5. For large values of K, regular settings approach continuous settings which deserve a 
special derivation. 
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r 


gtffr) : M 


gf (r) : M 


g^O) : M 


ef (r) : M 


g¥(r) : M 


n 
u 


1 
L 


COS ZCt 


tub ^±a 


PAC /I /~V 

COo 


cub za 


i 

_L 


cos 2a 


1 




pAo Ar\i 


pAc /I r\i 


o 
& 


/"> A S~\r 


lud za 


1 


lud za 


pac /I r\i 

cub ^ta 


Q 
O 


f~*{~\£! A r\i 


PAC /I y~V 
tUD 4:Ct 


COo ZCt 


1 


PAC / f~\i 

cob za 


4 


cos 2a 


r*/~iC /1 r\i 


phq A\rv 


pAQ / f\l 


1 

_L 


£ 
o 


i 


cud za 


PAC /I 


coo ^a 


cos 2a 


6 


cos 2a 


i 


cos 2a 


cos 4a 


cos 4a 


7 


cos 4a 


cos 2a 


1 


cos 2a 


cos 4a 


8 


cos 4a 


cos 4a 


cos 2a 


1 


cos 2a 


9 


cos 2a 


cos 4a 


cos 4a 


cos 2a 


1 



Table 10: Gauge distributions gf?(r)/M for 5 settings M = (1/2) sina ; a = tt/10 

2.3.5 Continuous 2D-EPR-B experiment 

For a continuous ensemble of settings, we have to construct a continuous ignition set (i.e., of 
cardinality Hi) and thus to define a convenient projection function. Next we will derive a continuous 
set of gauge probability distributions. 

Ignition set For any setting labelled 9 let (£|6>) be the local target £ € X, given 9 £ 0. For 
any real A, define the ignition state A by the intersection of local targets: 

A = p| {(U(6,X)\6)} 6 A, (28) 
ee[o,7r] 

where according to Eq. (|14[) . the binary digit £ 6 X is computed by a projection function 11(9, A). 
It is convenient to define the following dyadic square wave function of 9 and A guessed by analogy 
with Eq. (f25)) : 

11(0, A) = i[l +sgncos(0 — A)]. (29) 

Then, we will replace A by A when necessary and we will use a subset of A defined by Eq.(|28[) and 
corresponding to < A < 27t. This leads to a single solution. When compared with the regular 
discrete settings, we find a close similarity. This justifies the name of trigonometric order for the 
double-plateau functions in Sec. 12.2.11 

Gauge probability distributions For each 9 e & we are going to derive a gauge distribution 
density, gs(A)dA. Define the interval 

I(x,9) = {A j 11(0, A) = x} 

and let J(xo, xi, uo, m) = I(xo,uo)f]I(xi,ui). To compute the solution, let no = 9. For each 
setting ui £ © and each outcome pair xo, xi, we have the functional equation: 

Vui, Va;o,Va;i : / g e (A)dA = P(a; ; Ki|«o = 9; in) (30) 

where the unknown function is ge(A) or simply, ge(A). The solution has to be computed for each 
interval. For instance, for u± > 9, we have J(l, 1, 9,u\) — [ui — 7t/2, 9 + it/2] and P(l; 1\9; ui) = 
(1/4)[1 + cos(# — Mi)]. Finally, the result is the following working set of gauge distributions [25l |: 

g e (A) = (1/4)| cos(0-A)| (31) 

Eqs. (I29|l and (|31[) are closely related with Eq. (|25f) . Now the classical collapse is described by 
the process of Sec. 12.2.31 with some straightforward changes. Clearly, this model exactly emulates 
the 2D-conventional EPR-B experiment and leads to a violation of Bell's inequalities. 

3 Simulation of multipartite entangled systems 

We will now extend the above results beyond the particular case of Bell-type system to general 
systems with any number of regions. For simplicity, we will only deal with dichotomic outcomes. 
This means that we only consider an ensemble of qubits when the system mimics a quantum 
mechanical situation. In any case, the concept of local consistency remains the key property of 
classical analogues of quantum systems. This condition allows only nonsignaling correlation in the 
formulation of Popescu and Rohrlich [3^ |. 
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3.1 Classical analogues of quantum systems 



Consider an extended object in equilibrium composed of n correlated entities {£o, <£i, ■ • ■ , 
located in n space regions TZo, IZi,. . . , lZ n -\. The equilibrium can be broken and the future 
outcome Xi will depend randomly upon both the symmetry of the object and the free choice u, of 
the observer. When m = 6k, it is convenient to give a name to the pair (i, k) of one setting in one 
region: 

Definition (Configuration). We will call configuration a pair (i,k) of one region TZi and one 
setting 8k- Each configuration is uniquely labelled by an index 7 = k + iK , 7 £ [0, nK — 1], 

Let x be the outcome vector (xo; x\; . . . ; x„-i) and u the setting vector (uo;ui, . . . ; We 
will name again local target the entry (xi\m) and global target or simply target the pair (x|u). Let 
X = {(xju)} be the set of all targets. Depending upon the symmetry of the object, each target has 
a given conditional probability P(x|u). There are card(T) = 2™ x K n = (2K) n distinct targets, 
subjected to K n normalization constraints, namely: 

1 1 

••■ ^2 P(a;o;a;i;...;a;„_i|uo;wi;...;« n -i) = 1 

2:0=0 x T1 _i— 

When the system mimics a quantum situation, each setting-vector u can be regarded as an 
orthonormal basis of a Hilbert space H whose unit vectors are labelled by the 2 n outcome-vectors 
x. Therefore, each target (x|u) may be identified with a unit ket |x,u) of the Hilbert space. 
When the system is in the reduced state p, where p is a density operator, the target probability is 
computed as 

P(x, u) = Tr(p|x, u)(x, u|) = (x, u|p|x, u) 

More generally, given an extended object in equilibrium defined by a probability distribution 
P(x|u), we aim to describe classically a break of equilibrium of the system. Beforehand, we will 
clarify the concepts of local consistency in this context. 

3.1.1 Nonsignaling correlations, complete local consistency and degrees of en- 
tanglement 

We no more necessarily require total correlation, i.e., equal outcomes for identical settings in differ- 
ent regions, and accept any degree of entanglement. By contrast, we do conserve local consistency 
in order to allow independent partial measurements. This means that a definite local probability 
holds in each region, irrespective of the settings in the other regions. More generally, we will define 
the concept of complete local consistency when it is possible to derive consistent subsystems inn-r 
regions simply by ignoring r regions. Again, this condition only allows nonsignaling correlations 
between any pair of regions. A limit case concerns independent regions, when the probability of 
the whole system is simply the product of the probabilities of each region. Then, the system is 
separable and the entities are no more correlated. 

Definition (Complete local consistency). A set of n entangled classical entities {£0, £1, • ■ • , 
£ n _i}, is completely local-consistent when the measurement in one or several regions does not 
affect the marginal probabilities in the other regions. 

For any partition of {0, 1, . . . , n— 1} into two subsets, say (after reordering the indices if neces- 
sary), {0, 1, . . . , r— 1} and {r, r+1, . . . , n — 1} with 1 < r < n — 1, irrespective of u r , u r +i, . . . , 
we have 

1 1 

Prob(a;o;...;av-i|wo;---;Mr-i) = • • ■ ^ P(a;o; ■ ■ ■ ;a;„_i|«o; . . . ;w„_i). (32) 

x r - x n _i=0 

In words, Prob(a;o; • . . ; x r -i\uo; . . . ; it r -i) only depends on the local settings, («o; . . . ; u r -i). Thus, 
the free choices of u r , u r +i, ■ ■ ■ , u n —i and the subsequent measurements of x r ,...,x n -i do not 
affect the probabilities in regions IZo, TZ\, . . . , lZ r -i, as far as these choices and the subsequent out- 
comes are ignored. It will be convenient to name partial target the entry (xo; . . . ; a; r -i|tto; . . . ; u r -i)- 
Therefore, when complete local consistency holds, any partial target has a definite probability. 
Note that when the n-region system mimics a quantum situation, Eq. (|32p describes the partial 
state in the subsystem {IZo, . . . , TZ r -i}, obtained by partial tracing over the subsystem {lZ r , ■ ■ ■ , 1Z„-i}, 

Conversely, it is possible to increase the number of regions simply by adding successively 
new regions. This is similar to the process of quantum purification. Any subset of regions is 
consistent and the global and local viewpoints are compatible. Then, a multipartite system can 
be decomposed into its subsets without local perturbation. Note that computation of the partial 
probabilities Prob(a;o; ■ • ■ ; x T -\ |ito; ■ • ■ ; u r -\) implies that the regions lZ r , !Z r +i, ■ ■ ■ are ignored. 
As stressed in Sec. (|2.1.2[) , partial probabilities are subjective by definition. 
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3.1.2 Degrees of freedom 



Besides the K n normalization constraints, complete local consistency provides a number of addi- 
tional constraints between the target probabilities P(x|u). Since any subset of regions is consistent, 
it is possible to enumerate the degrees of freedom in each m-partite subsystem in isolation. In 
each single region, we have K local independent parameters, e.g. Prob(0|«i), because Prob(l|«i) 
is derived by normalization, and thus nK degrees of freedom for n regions and 2K for 2 regions. 
In each pair of regions, we have seen in Sec. 12.2.21 that we have K + 2K degrees of freedom. 
Since 2K local parameters are already counted we have K 2 joint bipartite degrees of freedom, e.g., 
Prob(0; 0\v,i\ Uj). As a result, we have {^)K 2 bipartite degrees of freedom for the general n-region 
system. Similarly, it is easily shown by induction that there are K m joint m-partite degrees of 
freedom for a general m-region system and hence {^)K m m-partite degrees of freedom for the 
general n-region system. Finally, in the general multipartite system, the number of degrees of 
freedom is 

Because only nK degrees of freedom are local, in general, (K + 1)" — nK — 1 degrees of freedom are 
shared. These results are summarized in Table fTTl Note that only the total number of constraints is 
meaningful, and for n > 1 one may prefer to consider that there is only 1 constraint of normalization 
and 2 n K n — (K + l) n constraints of local consistency. One or several degrees of freedom can be 
suppressed by additional constraints. For instance, in order to force all degrees of freedom to be 
shared, besides the constraints of normalization and total consistency, we have to add a number of 
nK new constraints. For n — 2, these 2K new constraints may be P(0; l\6k',6k) = P(l; Q\0k\ 8k) = 
forcing total correlation. Similarly, for n > 2, we can suppress both only-local and only-bipartite 
correlations with a number of nK + (") K 2 additional constraints. 



Number of regions 


n 


Number of settings 


K 


Number of targets (x|u) 


(2K) n 


Total number of local targets 


2nK 


Constraints of normalization 




Constraints of local consistency 


K n {2 n - 1) - (K+l) n + 1 


m-partite degrees of freedom 


( n )K m 


Total degrees of freedom 


(K + 1)" - 1 


Shared degrees of freedom 


(K + 1)" -Kn-1 


Number of gauge distributions g 7 (Aj) 


nK 


Rank of the gauge linear system 


2(K+l) n - 1 



Table 11: Degrees of freedom in general multipartite locally consistent systems. Gauge parameters refer to one-step 
collapse. 



Specially, when the system mimics a quantum set of n spins, the density operator p is defined 
by 2 ™ — 1 real coefficients. On the other hand, the number of degrees of freedom of P(x|u) is 
(K + l) n — 1 ( Table [TTj) . Therefore, a given reduced state p may be uniquely described by a locally 
consistent system of n regions and K different settings where K — 2 2 — 1 = 3. This number of 
independent settings, K — 3, is clearly related to the fact that spinors belong to a representation of 
the rotation group in R 3 . As a result, we may conjecture that a number of properties in quantum 
systems may be proved simply in a classical model with just three independent settings. 

However, in the present model, we will accept any number of independent settings. Thus, if 
a locally consistent system of n regions is defined by a set of independent target probabilities for 
K > 3 settings, and if the system has to mimic a quantum situation, the settings are subjected to 
a number of compatibility conditions (in addition with other constraints in form of inequalities, 
e.g., Eg. [531 below). 

The degrees of freedom are only relevant for general systems. Indeed, a particular system is 
completely definite and has thus no free parameter. It is useful to define instead a number of 
coefficients of entanglement, characterizing the behaviour of the system subject to splitting. The 
main issue, arguably surprising, is the emergence of randomness. 



3.2 Randomness and entropy 

We regard the classical equivalent of a stationary quantum system as a spatially extended object 
in equilibrium. As long as the dynamic equilibrium holds, the system evolution is completely 



15 



deterministic. The potential probabilities P(x|u) describe only the symmetry of the object, exactly 
as a spinning roulette can be characterized by the prior probabilities of its wheel pockets, but before 
any break of equilibrium, by definition, the entropy remains zero. Similarly, in quantum physics 
the von Neumann entropy of a pure state in unitary evolution is zer 

Consider for a start the case of bipartite systems. Randomness can arise from two causes, 
measurement and partition of the system. 



3.2.1 Entropy of measurement in bipartite systems 

Firstly, randomness arises by measurement, because the initial equilibrium is broken. This is a 
well known situation in classical physics: After a break of symmetry, several final outcomes are 
possible at random. The randomness is then described by the conventional Shannon entropy 
accounting for all constraints at the very best. Since the break of equilibrium follows a free choice 
of the settings, there are as many entropy coefficients as possibilities. For a finite number of K 
settings in the 2-region context, the entropy is described by a square isT-matrix, Si(«o,ui). 



Global measurement For a global measurement, we have, 

1 1 

Si(uo,iti) = 2Z ^2 — P(xo; xi\u ; m) log 2 P(x ; Xi|uo; mi) (33) 

xq=0 xi=0 

We will name Si global measurement entropy. The point 'global' is a matter of context. 

Note that this function Si characterizes the classical randomness following a conventional break 
of equilibrium. When the system mimics a quantum situation in a reduced state p, Si should not 
be mistaken for the von Neumann entropy, — Tr(plog p), which characterizes the 'purity' of the 
state and is not related with its the amount of informatiorQ 



Partial measurement It is possible also to consider only a partial a measurement, e.g., in 
region TZo, and we have 

1 

Si(ito) = 2Z -Prob(xo|«o) log 2 Prob(a;o|wo) (34) 

XQ=0 

where Prob(a;o|Mo) is computed from Eg. ((71 13"2)) or Eq. (|21|> . We will name Si partial measurement 
entropy. Due to local consistency, the point 'partial' is again a matter of context. 

When the system mimics a quantum situation in a pure state p, this function Si is closely 
related to the so-called 'entanglement entropy' E\ in quantum information [45l |. namely e.g., in 
region TZ a , 

Ex{no) = -Tr (p logpo) (35) 

where Trj(.) stands for partial tracing over the region TZi and pi — Tri_i(p), (i 6 {0, 1}). Still for 
pure states, we have [46l |. 

Ex{K ) = E 1 {'R, 1 ). (36) 

In the present model, owing to Eg. 1101 the same equality Si(uo) = Si(ui) applies to Bell's-type 
systems when Uo = u\. However the equality does not hold in general. For separable systems, 
E\(1Zq) and Ei(lZ\), or Si(ito) and Si(ui) in the present model, are even independent. What is 
most important, neither the function E\ in Eo. (|35|l nor Si in Eo. (|34ll do grasp the level of entangle- 
ment. Therefore, we regard the partial measurement entropy, Eo. (l34|) . as completely distinct from 
the entanglement entropy, even if the two values coincide for bipartite totally correlated systems, 
and we will use a different definition for the 'entanglement entropy' in the next section. 

6 Actually, the von Neumann entropy is equal to zero when the reduced state of the system is pure, i.e., not mixed 
(See Sec l3.2.1T l. This implies that the evolution is deterministic. 

7 In the present model, if we are concerned by the degree of 'purity', i.e., the correlation of the outcomes for identical 
settings in the two regions, a witness function could be 

1 

c(6) = ^2 -P(x; x\0; 9) log 2 P(x; x\0; 9) 

x=0 

With the current convention (see footenote [TJ we have c(6) = for pure states and this function is closely related to 
the von Neumann entropy, — Tr(plogp). With the opposite convention, used in Sec. [2] we have c(9) = 1 bit for totally 
correlated systems. 
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3.2.2 Entropy of entanglement in bipartite systems 

Less expected, a second source of randomness arises in entangled systems without break of equi- 
librium. Indeed, it is possible to split the whole system into its two parts and consider each region 
separately. Then two cases are possible: either the two subsystems are independent, or the two 
regions are entangled. In the first case, the two regions are separable, all information is strictly 
local and no randomness appears. The entropy of the whole system is the sum of the entropy of 
each subsystem, i.e, the entropy is extensive [471 ] . 

By contrast, in the second case, local information is insufficient to completely describes the 
whole system and additional information is needed. By splitting, we decide to ignore all non local 
information and our description of the system becomes incomplete. In other words, the system is 
no more completely definite. Therefore, without any break of equilibrium, the entropy of the two 
separated subsystems is now positive, clearly greater than the zero entropy of the full system. As 
a result, entropy is no more extensive. 

In quantum theory, this is generally viewed as an astonishing aspect of entanglement. Actually, 
this situation is well known in classical thermodynamics in the face of long range interactions [3ll . 
|33 |. In the present case, the only possible paradox lies perhaps in the fact that the splitting is 
nonsignaling, and thus each region remains locally consistent. 

Quantitatively, we aim to evaluate non-local correlations between a pair of local targets, (xo\uo) 
and (xi\ui), where Xi is regarded as a random variable while Ui is viewed as a parameter. A similar 
issue holds in the conventional Shannon transmission theory. The concept of mutual information 
was indeed devised to evaluate the correlation between the data at the two ends of a communication 
link. This concept was extended to multipartite systems by R. Fano [48J. In this background, 
the entropy of measurement Si can also be called self information. For ease of exposition, it is 
convenient to use a shorthand Xi for the local target, (xAui). Then (Xq; Xi) describes the global 
target (xo;xi\uo;ui). Mutual information is defined as [49||, 

I(X ;X 1 )^ ) E[lo g2 Pmh ^ X ^ i 



Prob(X )Prob(Xi) J 

where E[.] stands for the expectation value with respect to Prob(Ao;Xi) = P(xo',xi\uo;ui). Mu- 
tual information does not depend on the order of the two targets, 



Define, 



We compute easily, 



I(X ;X 1 ) = I(X 1 ;X ) (37) 



S2{uo,ui) = I[(xo\uo); (a;i|«i)] = I[(x!\ui); (x \uo)] 



s 2 (u , ui ) = Y, E P(*o;*i|«o;ui)iog a ( ^r^ferrn ( 38 ) 

' — ■ - — • \ Prob(a;o|wo)Prob(a;i|ui) ' 



Xq — X\= 



In the present context, it is convenient to call this coefficient S2 bipartite entanglement entropy. 
Depending upon the background, mutual information is also named (apart from the sign) relative 
entropy [501 1 or Kullback-Leibler divergence [5l|. 

Now, rename the local probability distributions Prob(sojiio) = Qo(a;o|uo) and Prob(a;i|iii) = 
Qi(xi|mi) and define a global distribution Q(x|u), as 

Q(x|u) <d = f) Qo(a;o|«o) x Qi(xi|«i) (39) 

Q = Q(xju) describes a separable system, derived from the original system P = P(x|u) by removing 
all non-local correlation. Using the relative entropy S(P||Q), 

1 1 

s(piiq) =EE p 0"° 

;a:i|uo, ui) logP(:ro; xi\uo, ui) - P(a;o; X\\ua, m) log Q(a; ; xi\u , «i) 

xq=0 xi=0 

we have identically 

S 2 (u) ( = f) S(P||Q). (40) 

In words, 82(11), is the relative entropy of the entangled distribution P(xju) with respect to the 
separable distribution Q(xju), Eq, (|39p . Therefore, S 2 (u) is always positive and is zero only for 
a pair of separable regions. It is invariant by interchange of the regions but not, in general, by 
permutation of the settings. Note that the relative entropy is non symmetrical by interchange 
of the relevant distributions, namely, the joint distribution P(x|u) and the separable distribution 
Q(x|u). 
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Let Qq (:ro | wo) and Q'i(xij-ui) be two arbitrary local distributions in IZo and IZi respectively. 
The distribution Q', 

Q'(x|u) ( = f> Qo(a;o|Mo) x Qi(xi|ui) 
is separable. Accounting for the definition of Qo and Qi we obtain, 

S(P||Q) - S(P||Q') = -S(QollQ'o) - S(Qi||Qi) < 0. 

Therefore, 

S 2 (u) = S(P||Q) = min S(P||Q'), 
Q'er> 

where T> is the set of all separable distributions. The minimum is clearly obtained for Q' = Q. 

Furthermore, if the two regions are totally correlated, the whole information is shared be- 
tween the two parties and the entanglement entropy for a given setting, S2, is equal to the local 
measurement entropy, Siioc, of any region. 

When the system mimics a quantum situation in a reduced state p, the separable distribution 
Q mimics a separable system defined by a density operator a, 

a — (Jo <8> err with cro = Tri(p) and o\ = Tro(p), (41) 

where the partial trace Tr;(.) stands for Tr-]z i (.) (i € {0, 1}). The quantum relative entropy of p 
with respect to a is [52| . 

S(p||a) = Tr(plogp)-Tr(ploga), 
while the quantum entanglement entropy, -E 2 (p) was defined by Vedral et al [HJ as, 

£ 2 (p) ( = f) minS(p||a'), (42) 

<t'€E> 

where T> is the set of all disentangleojf^tates. Clearly, if the system is separable, p — a and E2 = 
and conversely. As a result, E2 — if and only if the system is disentangled. More generally, with 
the same proof as in the classical case, it can be shown that minS(p||<r') is obtained for a' — a. 
Then -E 2 (p) — S(pjjcr). Therefore, in bipartite systems, the present definition, Eq. (|38l 140 [I . is 
closely related to the definitiorJf|of the quantum entanglement entropy Eq. (1 12[| [54H-56l |. We will 
revisit this concept for n-partite systems. 

For pure state, S(p||<r) coincide with the von Neumann entropy of the partial density pi in any 
region TZi. 

S(p||o") = -Tr (po logpo) = -Tri(pilogpi), 

but in general, this relation does not hold. 

In the present model, for a discrete classical system with K possible settings, the entanglement 
entropy S 2 (u) is a square matrix of order K, in general non symmetrical. This result is not 
surprising because we know that we have K 2 bipartite degrees of freedom (Tables [1] and lll[) . 
Therefore we can identify these K 2 parameters with the K 2 coefficients of bipartite entanglement. 
As a result, a bipartite system is separable if and only if its entanglement entropy matrix is zero. 

The bipartite entanglement entropy S2 measures the amount of information required to account 
for non-local bipartite correlation. For instance, in a previous paper [26], we have described a dice 
game between two parties, Alice and Bob, to illustrate the concept of classical entanglement. In 
this example, it is easy to compute an improbable 'entanglement entropy' between the two players 
as a 3-matrix, given in Table [12] 

For a continuous ensemble of settings, the entanglement entropy is a function 82(1*0, Wi). Spe- 
cially, for a pair of EPR totally entangled entities, the entanglement entropy with respect to a pair 
of settings 6 a and 9b is 

S 2 (9a,9 b ) = (1/2) [(1 + cos6»)log 2 (l + cos<9) + (1 - costf) log 2 (l - cos<9)] 

where 8 — 8 a — #(,. Clearly, S2 < 1 and the maximum is equal to 1 bit, obtained for 6 — 0. This is 
just the amount of information required to transmit classically the non-local correlation between 
the two entities of the pair [5?| . 

8 We take always the word 'disentangled' as synonymous of 'separable', while in reference |5''I the signification is 
complex and more subtle. However, the difference seems not very clear. 

9 The definition Eq. 11381 is constructive and provide a necessary and sufficient condition of classical bipartite entan- 
glement. We can obtain similarly a constructive definition of quantum bipartite entanglement by droping 'min' and 
substituting a for cr' in Eq. ll42t . 
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S 2 (A a ,A b ) 


Ai 


A 2 


A 3 


Ai 


1 


0.35 





A 2 


0.35 


1 


0.35 


A 3 





0.35 


1 



Table 12: Example of classical entanglement entropy (in bits) between two players, Alice and Bob, in a totally correlated 
classical dice game described in Ref. |2(j|. The game is implemented with three dice labelled Aj, [k = 1,2 and 3). At 
each run, each player selects freely one die, A a and Aj, respectively. The coefficients measure the amount of information 
required to account for non-local correlation. 

3.2.3 Multipartite entanglements 

When considering multipartite measurements, there is no new topic. By contrast, splitting of 
the system raises different questions because there is a number of different ways to group the 
regions together. Again, the same issue is encountered in communication networks: In transmis- 
sion theory, the natural extension of the bipartite mutual information is the multivariate mutual 
information [48| . Similarly, we will define multivariate mutual information coefficients between 
local targets, computed recursively as [4^|, 

I(Xo; Xi ; . . . ; X n -\) ( = ' I(Xo; X\ ; . . . ; X n -s; X n -2) 

- I[{X ; Xv, . . . ;Xn-s;X„- 3 )|X„-i)] (43) 

In Eq. (|43p . Xi is still a local target, shorthand for (xi\ui), where Xi is a random variable and m a 
parameter. It can be shown that multivariate mutual information is completely symmetrical with 
respect to Xo,Xi, . . . ,X n -i even if this is not manifest in the definition. Contrary to bivariate 
mutual information which is always positive, multivariate mutual information can be either positive 
or negative. 

Shannon's information theory has a nice set-theoretic structure 49, 58, 59] recalled in Table [T51 
Thanks to this structure, computation of entropic expressions is simplified by a graphical repre- 
sentation, named information diagram and similar to the conventional Venn diagram. Entropy or 
mutual information / are associated with a signed measure fi. For simplicity, and without loss of 
generality, we will use the same symbol / to describe both entropy and information (while entropy 
is usually symbolized by H). Multivariate mutual information, Eq. (|43[) , is directly associated with 



Information theory 


Set theory 


I(X) 


H(X) 


(Xv,X 2 ) 


X\ n x 2 


(X U X 2 ) 


X\ u x 2 


(Xi|X 2 ) 


X\ — x 2 



Table 13: Set-theoretic structure of Shannon's information theory. Each variable X is represented by a set X. Entropy 
or mutual information / are associated with a signed measure fi. Information expressions, like I(X\, X2), 1{X\ : .X'2) or 
I (X\ X2) are respectively associated with fJ.(Xi U X2, n(Xi n X2 and fi(Xi — X2) (where X\ — X2 = X\ n X|). 

I(X ;...;X n - 1 ) = l x(X n---nX n - 1 ) (44) 

and thus I(Xo; . . . ; X n -i) is clearly completely symmetrical with respect to Xq, . . . , X n -i. For 
instance, when n — 3, Eq. (|43|l is represented by the information diagram shown in FigfJ] Given 
a setting vector u, each local distribution Xi is represented by an area and I(Xi) is the partial 
measurement entropy, or the self-information of this variable. We consider only binary outcomes 
and thus I{Xi) < 1 bit. For instance, Eq. (|34|) . describes the self-information I(Xq). The n 
regions delimit a number of 2 n — 1 atoms of information, fj,i, . . . , fj,2 n -i- For n = 3, we have 
I(X ;X 1 ;X 2 )=fJ, 7 . 

In practice, it is convenient to derive first the measurement entropy I(Xi a , . . . , Xi r l ) of each 
partial subsystem (lZi , ■ ■ ■ ,TZi r ) for r = 1, ...,n, because each probability distribution is easily 
computed. For instance, if {io, . . . , i, — 1} = {0, . . . , r — 1}, the relevant partial probability in 
the r-region subsystem is given by Eq. (|32p . We have 2 n — 1 such partial subsystems. Next, 
we express the measure /ii of each atom of the information diagram (Fig. [T]) in terms of these 
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measurement entropies. We have 2 n — 1 atoms and therefore we obtain the mutual information 
I(Xq;Xv, . . . ; Xn—i) = /ti2™-i by solving a linear system of 2™ — 1 equations for 2™ — 1 unknowns. 
For example, when n = 3, we derive easily the following equation by inspection of Fig. [TJ 
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and we obtain 
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The last row reads 

(i 7 = I(X ; Xi;X a ) = I(X ) + /(-Xi) + I(X a ) 

-I(X ,Xi)-I (X , X a ) - I(Xi , X„) + I{X , Xi , X a ) (45) 

For instance, it is possible to mimic the 3-region GHZ-state with the setting (Z; Z; Z), where Z 
is the axis of common entanglement. Similarly, we can mimic the W-state with the same setting. 
The measures (n of all atoms are given in Table 1141 
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0.667 





0.667 


0.667 


-0.415 


1.170 



Table 14: Signed measures fii and total entanglement S = 2(17 + /13 + ^5 + /i6 (in bits) for the tripartite GHZ and W 
states with the setting (Z; Z; Z). 

More generally we have, 

I(X ; Xi; . . . ; X n _i) = £ I(X t ) - g I(X t , X } ) + 

J2 I(Xt,X 3 ,X k ) -■■■ + (lr-'liXo, . . . ,Xn_i) (46) 

i,D,k 

where the summation is taken over all possible combinations of subscripts and, 
1 



I(Xi) = Si(Mi) = V> — Prob(a^ijui) log 2 Prob(xi|«i) 



^ -P(x,u)log 2 Prob(xj|Mi 

xq=0 x T1 _i— 



x j —0 X.: — 



I(Xi, Xj) = Si(m,Uj) = ^2 — PTob(xi\Xj\ui\Uj)log 2 Prob(xi;Xj\ui-,Uj) 

1 1 

= yZ ■■■ 5Z -P(x,u)log 2 Prob(»j;a;j|ui;Wj 



xq— x n _i— 

etc. In the present context, we will use the following terminology: 
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Definition (n-partite entanglement entropy, degree and coefficients of entanglement). We will 
name n-partite entanglement entropy the expression 

Sn(u) = J(X ;Xi;...;jr n -i). (47) 

When different from zero, S„(u) defines one degree of n-partite entanglement, e n = 1. We will 
call entanglement coefficients its K n components. 

The K n entanglement coefficients can be identified with the K n n-partite degrees of freedom 
of a general n-region system. As a result, if M n 7^ 0, we will regard the n-region system as 
n-partite entangled. In general, S n (u) is not invariant by permutation of the settings but the 
object describes the n-partite entanglement as a whole and does not depends on the order of the 
n regions. For n = 3, we have 



ss(u) = E E E p w«) 



xq— x\— X2— 

/ ¥roh{x Q \x\\uo\U\)¥Tob(x 1 \X2\u 1 \U2)¥roh(x2-,x \u2]u ) \ 
2 \P(io; x\\x2\uo; ui; U2)Prob(a;o|wo)Prob(a;i|ui)Prob(a;2|«2) / 

For instance, from Table [TH we have S 3 (Z; Z) = 1 bit for the GHZ state and S 3 (Z; Z; Z) = 
-0.415 bit for the W-state. 



/'2 



M3 / \ A*6 



Mi 



/ M7 \ 












\ M5 / 







Figure 1: Information diagram of three parties, Xo, and X2. Using Table [T3l we have for instance I(Xq) 
Mi +M3 +M7, i"(X ;Xi) = At 3 +/t7, /(X ;Xi;X 2 ) = 117- 



3.2.4 Characterizing entanglement 

Definition (Degree of entanglement). When a n-region system is n-partite entangled, we will say 
that its degree of entanglement is n and set e n =l. Otherwise, its degree of entanglement is the 
larger number of regions, m, which constitutes a m-partite entangled subsystem. 

Finally, any given system has a particular structure described by the scheme (e.2, e%, ■ ■ ■ , e n ) of 
its bipartite, tripartite,. . . , n-partite degrees of entanglement. The number of m-partite coefficients 
of entanglement non identically vanishing is thus e m K m , with e m < ( n ). 

Definition (Entanglement scheme). Given a particular n-partite system with K settings, we will 
call entanglement scheme the sequence (fii,ez, . . . ,e n ) of its bipartite, tripartite,. . . , n-partite de- 
grees of entanglements, where e m < f n ) . The overall number of entanglement coefficients is 



N = E e™K r - 



m=2 



In addition, we will define later K n total entanglement coefficients by Eq. (|50H in Sec. (|3.2.5jl . 

Given the probability distribution P(x|u) the coefficients of entanglement are easily com- 
puted. For instance, Bell-type systems are bipartite, totally correlated systems, characterized 
by the scheme e% = 1 (with i^)K 2 = K 2 coefficients of entanglement). Among tripartite sys- 
tems (Sec. 1431 below) . GHZ-states (considered with only two settings) are characterized by the 
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scheme (e.2, e%) = (0,1) and K — 2, with (g)7f 3 — 8 coefficients of entanglement, while general 
W-states (with two settings) are characterized by the scheme (62,63) = (3,1) and K = 2, with 
(2) A" 2 + K 3 = 20 coefficients of entanglement, but only 2 tripartite coefficients (up to 8) are non 
zero. Separable systems are trivially characterized by = e$ = • ■ • = e n = and thus zero 
coefficient of entanglement. 

The current results are rather paradoxical. We have elaborated upon a close analogy with 
signaling networks except that the present system is . . . nonsignaling! However that may be, 
information theory has established a number of helpful inequalities between entropy functions, 
which are of course currently valid. 

3.2.5 Information inequalities 

In quantum information theory, entanglement in multipartite systems is not yet fully understood. 
By contrast, in the present model, the information of all atoms is easily computed and we have a 
number of tools at disposal stemmed from information theory. We recall first the polymatroidal 
axioms. Next we define the notions of 'total entanglement' and maximally entangled systems. 

Polymatroidal axioms Let a, ft be two subsets of [0, n — 1]. If a = {ii, 12, . . . , i r }, define 

7(a) ( = f) I(X n ,...,X lr ) 

(for simplicity, we conserve the same symbol 7). 

The following polymatroidal axioms [SsJ describe the Shannon information inequalities. 

1. 7(0) = 

2. If a C P then 7(a) < 1(0) 

3. 7(a) + I(/3) > I (a n 0) + 7(a U 0) 

This tool is powerful to derive important entropic relations. Note that this is not the complete story 
since there are inequalities not derivable from this set of axioms, called non-Shannon inequalities. 
We will just give one example, and define the concept of total entanglement. 

Total entanglement In transmission theory, Watanabe [dQ| | has defined a global coefficient 
£, named total correlation as, 

n-l 

£(X ; Xi', . . . ; X n -i) = Y, J ( X - H^Xx,. . . ,X n _i) 

i=0 

For instance, in tripartite systems (Fig. [l]), we have 

£ (X - Xr, X 2 ) = 2/i7 + fj, 3 + fj, 5 + fi 6 

It is possible to express £(Xq; Xi; . . . ; X„_i) in terms of relative entropy as follows: Given a setting 
vector u, construct a separable probability distribution Q = Q(x|u), as 

n-l 

Q(x|u) < = ) J]Prob(x 4 |«i) (49) 

i=0 

Now we have, with P = P(x|u), 

£(u) = S(P||Q) 

As a result, total correlation is always positive and is zero if and only if all regions are independent. 
In the present context, it is relevant to use a special term, 

Definition (Total entanglement). Given a particular n-partite system with K settings, we will 
call total entanglement the function 

£(u) ( = f) S[P(x|u)||Q(x|u)]. (50) 

Note that the concept of total entanglement can be easily translated into quantum formalism. 
Consider a quantum n-region system in the reduced state p and let Oi be the reduced state in 
region IZi, obtained by tracing over the n — 1 other regions. Let a be the the separable state 

n-l 

Define the 'quantum total entanglement' as 
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£(p) = S(p\\a) (51) 

Clearly, when all regions are separable, £(p) — . We have already noted (Eq. 142 f) that the 
quantum entanglement entropy, E(p) was defined by Vedral et al as, 



E{p) <d = f) rninS(p||a'). 

o'€T> 

where T> is the set of all disentangled states. It can be shown again that the minimum is obtained 
for a' = a, then E(p) — S(p\\a) = £{p). As for bipartite systems (footnote [5] above), this provides 
a necessary and sufficient witness of entanglement. 

In the present model, since we only consider binary outcomes, we have 

I(Xi) < 1 bit 

As a result, we have the inequality < £ (u) < n bits. Actually this maximum can be improved to 
n — 1 bits. The proof is trivial for n = 1. Now, thanks to the polymatroidal axioms, the addition 
of one region X m can only increase the total entanglement by I(X m ) < 1 bit. Therefore, we have 
for a n-region system, 

< £(u) < n- 1 bits. (52) 

The upper bound can be saturated, e.g., when the measures pi of all atoms of information are 
zero, except for I(Xo; Xi; . . . ; X n -i) = 1 bit. This correspond to the general n-region GHZ-states, 
when the same setting, corresponding to the axis of entanglement, is used in each region. 

Definition (Maximally entangled system). We will say that a n-region system is maximally en- 
tangled when its total entanglement for at least one setting-vector is equal to n — 1 bits. 

For instance, among bipartite systems, EPR-pairs and specially the singlet state are maximally 
entangled (Table [17b below) as well as the so-called PR-box (Table [18b below), with £ = 1 bit. 
Among tripartite systems, the GHZ-state with the setting (Z; Z; Z) is maximally entangled with 
£{Z; Z; Z) = 2 bits (Table [Till. Note that with only two settings X and Y, the GHZ-state is not 
maximally entangled (Table [2Tb below). The W-state is not maximally entangled with a maximum 
entanglement of 1.170 < 2 bits (Tables [HI and [T9"b below). 

This framework is likely to be helpful to discuss a number of issues like 'entanglement monogamy' [6 
or collapse mecanism (Sec. 13.3. II below). These discussions are beyond the scope of this paper. 

To sum up, we have a necessary and sufficient condition of entanglement, namely, that total 
entanglement £(u), (Eq. [50}, is non zero for at least one setting-vector u and a bound of maximal 
entanglement, namely, £(u) < n — 1 bits for a n-region system. 



Classification of locally consistent systems In quantum information theory, an impor- 
tant challenge is to characterize entangled and separable systems. A number of entanglement 
witnesses are known. Bell-CHSH inequalities were first defined for only two regions. In 1990, 
David Mermin [62T | constructed a set of multipartite inequalities similar to the CHSH expressions. 
Since then, a large range of witness operators have been described [HE! ]. In the present model, all 
these operators are valid, but the concept of total entanglement provides an easy necessary and 
sufficient condition of entanglement. 

Beyond the separation between entangled and separable systems, a second dividing line may 
be drawn between quantum and 'super-quantum' systems, i.e., systems that can or not mimic a 
quantum mechanical situation [36]. According to Tsirelson, the characterization of quantum state 
was first raised by Anatoly Vershik, but the very question was addressed quantitatively by Boris 
Tsirelson himself [63|, |64j, who showed that the CHSH inequality, Eq.JBJ, became for quantum 
systems, 

|S(A, B) + S(A', B) + S(A, B') - S(A', B')\ < 2V2. (53) 

In addition, Tsirelson identified this bound, 2\/2, as a special value for two regions of the Grothendieck 
const anlF°l. defined in general topological tensor product spaces 65]. In terms of Hamming diver- 
gence, Eq.{T}, we have also, 

1 - V2 < d(A, B) + d(A', B) + d(A, B') - d(A', B') < 1 + V2 

The generalization of Eq. ([(J I53[) to multipartite systems was addressed by Mermin [r33 | but 
the surprising result was that the same bound is valid for both quantum and 'super quantum' 
systems. We will give in Sec. 13.3.21 a criterion based on the evaluation of the Tsirelson bound 
following a process of classical multipartite collapse. 

10 The algebraic nature of the Grothendiek constant is unknown and so, the constant is generally qualified of 'enigmatic', 
possibly because Alexandre Grothendiek is viewed as an enigmatic mathematician. 
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3.3 Classical multipartite collapse 



Consider a n-region classical system in equilibrium. The equilibrium can be broken as follows: In 
one or several regions IZi, an observer, Oi, (1) selects freely a setting Ui, element of a given set © 
of K different settings, (2) breaks the local equilibrium in a way reliant on u;, and (3) measures 
in IZi an observable Xi, element of a set X = {0, 1}. 

3.3.1 Collapse mechanism 

We need first to discuss the 'collapse mechanism', in the sense of 'reaction mechanism' in chemical 
kinetics. Indeed, the global process defined for Bell-type systems turns out to be insufficient to 
describe all multipartite collapses, even for independent regions. 

The EPR paradox was cleared up by a direct gauge implementation of the probability distribu- 
tions in the two regions. Actually, a gauge selection defined a leading region, say IZo, but afterward 
the two regions entered into the model on just the same footing. We have computed a probability 
distribution, </k (Aj), located on the boundary of the two regions, accounting for the free choice 
uo = ko in the leading region and governing both (a;o|fco) and (xi\ko;ui). 

However, local realism does not require that the same distribution governs the two processes. 
One can construct an indirect collapse as follows. (1) The settings (fco,fa) are sent to an ignition 
point X. (2) A gauge selection defines the leading region, e.g., IZo- (3) At the boundary of the 
two regions, the outcome xo = £o of the leading region is drawn at once, using Prob(a;o|fco)- (4) A 
gauge distribution Prob(a;i|£o; fcp; ui) is used in the second region IZi. Of course, thanks to the 
well known product rule [66j|, accounting for local consistency, the two processes give exactly the 
same final result, 

P(£o;:ri|fco;wi) = Prob(£ |fco)Prob(:n|£o; ko; ui) 

Clearly, the indirect process requires two successive trials. As a result, the number of gauge 
distributions is doubled. On the other hand, each distribution Prob(xi|£o; ko; u\) is simpler, and 
even trivial in bipartite system, because we have only to account for compatibility conditions in 
one region instead of two. Therefore, in the computation of the gauge distributions, the direct 
route can fail, while the two-step process turns out to be successful. As in chemical kinetics, this 
may be physically meaningful because the two-step route may require a delay, while the direct 
process can be viewed as quasi- instantaneous at the ignition point. 

In multipartite systems, we will proceed similarly. Thanks again to the product rule and local 
consistency, it is possible to define multi-step collapses, for instance, 

P(a;o;a;i;2;2;a;3|uo;wi;«2;W3) = Prob(xo|ito) x Prob(a;i|a;o;«o;wi) 

x Prob(x2|a;o;a;i;wo;wi;w2) x Prob(:r3|a;o; x%; xi\ u ; Mi; U2; u 3 ) (54) 

At each step, we have to select a new gauge region. For m steps, we have m gauge trials per run. 

Suppose, for instance, that several observers decide to proceed to a measurement and select 
freely a setting 0fc 4 . They send their setting towards an ignition point, located at the boundary of 
the relevant regions. A first gauge selection at the ignition point defines a leading configuration, 
say 70 = fco + Kio, i.e., a gauge region 7Zi and a gauge setting 6k - The first outcome £o is drawn 
using the local probability Prob(ajj |fco)- Now, if e.g., io = 0, the probability distribution in the 
remaining n — I regions is a gauge function 

Prob(xi; . . . ; x n -i\£a; feo;«i; • ■ • ; u n -i) = M x P(£o; Xi;... ;a;„_i|fc ; uv, . . . ;u n -i) (55) 

where M = M(£o> ko, v>i, . . . , u n -i) is a normalization factor. Note that Eq.(f5S]) only describes 
the updated probability of (x\, . . . ,x n -i) for the observer Co, because he has got alone a new 
information, namely uq = fco and Xo = £o- By contrast, the n — 1 other observers are not aware of 
what happens in region IZo, and their view of the n — I regions is still given by Eq. (|32|) . 

i 

Prob(a;i; . . . ; x n -i\uv, . . . ;u n -i) = ^ P(a:o; • ■ • ; x n -i\uo;uv, ...; u»_i). (56) 

xo=0 

Due to local consistency, the right hand-side of Eq. (|56[) does not depend on uo and we can plug 
fco in no. In general, Eq. (|55p and Eq. (|56p are different, because they describe different level of 
knowledge. 

At least formally, this process can describe a real quantum collapse. Therefore, if the original 
system mimics a quantum situation, the probability distribution Eo. (|55|) will also mimic a quantum 
system. More generally, it can be seen that Eq. (|55p is still locally consistent. Therefore we can 
define a second-step probability distribution P' as 

P'(a;i; . . . ;x n -i|i*i; • ■ • ; «n-i) ( = f Prob(a;i; . . . ;a; n _i|£o; k ; u\; . . . ; u n -i)- (57) 
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Owing to Eq,(154[>. this distribution P' describes now the n — 1 remaining regions (instead of 
P(x|u)). It is possible to iterate when wanted. Let m be the number of iterations. If m — n, the 
overall collapse is completed. Each observer gets his own outcome, Xi = £j but ignores the other 
outcomes. If m < n, the last gauge trial has to be performed in a subsystem of n — m regions. For 
this purpose, we will construct in the next section a convenient stochastic gauge system on the 
model of the Bell-type systems. If we fail to obtain a set of positive gauge probability distributions, 
this means that there is no locally realist collapse in m steps and we have to try with m + 1 steps. 
The last stage in just one region is trivial, and therefore, it is always possible to find a locally 
realist route whatever the number of regions and the degree of entanglement. In other words, 
given any locally consistent system, the global collapse can always be implemented by an arbitrary 
chain of n one-region collapses. 

Definition (m-step collapse). We will call m-step collapse a collapse composed of m-1 cascaded 
one-region trials and one final trial into a subsystem of n — m + 1 regions. 

In bipartite systems, we have always found by computer simulation that the direct collapse in 
just one step is feasible. This is not always the case in multipartite systems. We will see that a 
one-step process fails is some tripartite 'super quantum' systems. When more than three regions 
are involved, the one-step collapse may fail even when all regions are separable. 

Note that a n-step collapse, requiring a minimum of computational resources, can be helpful 
to simulate quantum algorithms. 

3.3.2 Using the Tsirelson criterion for more than two regions 

We have seen that any locally consistent system of n regions can always collapse by an arbitrary 
chain of n one-region collapses. Since a quantum system can formally follow this process, each 
step in the chain defines a quantum gauge subsystem of m regions (with m 6 [l,n]). Therefore, 
only a 'super quantum' system can lead to a 'super quantum' gauge subsystem. In other words, 
if one gauge subsystem is 'super quantum' we can conclude with certainty that the global system 
was 'super quantum'. 

Specially, evaluation of the Tsirelson bound allows the 2-region 'super quantum' subsystems 
to be detected. Beyond, multi-step collapses provide a criterion of non-quantum behaviour for 
more than two regions. For instance, we will use this criterion in Sec. 14.3.31 below to characterize 
a 3-region super-quantum system. 

3.3.3 One-step collapse 

In order to describe a classical one-step collapse in multipartite systems, we have to generalize the 
concepts defined in Bell-type systems, namely, ignition states, gauge probability distributions, and 
measurement process. 

Ignition states The ignition states will define the potential outcomes in the n regions. In 
Sec (|2.2.I|) we have constructed a fine-grained description of Bell-type systems. Now, we aim to 
construct a fine-grained description of multipartite systems by the same token. 

We have defined the local target (xi\v,i) = (£|$k)» m a region TZi by Xi = £ given Ui = 9k- 
When the system is completely local-consistent, the probability Prob(£|#k)i in a single region is 
uniquely defined. The number of distinct local targets is thus 2K for one region and 2nK for the 
whole. Next, we construct the subset Zms.y of X composed of all global targets compatible with 
one particular local target 

%i« fc )i = {( x l u ) / x i = £ ; ui = e k } 

Any local target (£|#fc)j in region TZi can be regarded as an union of global targets 

(£l**)i= U « x l u »- 

(x|u)6T (? | 9fc)i 

Conversely, any global target (x|u) can be viewed as the intersection of n local targets in regions 

TZq, . . . , TZn-l 

(x|u) = (a;o|«o) n {x\\ui) n • ■ • n (x n -t\u n -x) 

Let 3 = {j} be a set of integers. For i £ [0, n— 1], define n applications Hi(k, j) of [0, K— 1] x 3—)- 
X = {0,1} 

€i = ni(*,i). (58) 

We will call projection function in region TZi the function Hi(k, j). More concisely, it is possible to 
define a global projection function 11(7, j) m terms of configuration as 

n( 7 , j) ( = f) Ui(k,j) with 7 = k + iK (59) 



25 



Now, we construct a sequence of card(J) fine grains Xj by the following intersection of local 
targets: 

3 G 3 >-> Aj 

n-l A'-l 

h ( = f) n n (ni(fc,j)i^).. (60) 

i=0 fe i= 

There are if™ different setting vectors, labelled (ko, ki, ■ ■ ■ , fen-i)- Therefore, the maximum num- 
ber of distinct fine grains is 2 K . On the other hand, we have to define the outcome Xi for each 
configuration (i, fc) given j. There are only nK configurations 7 = k + iK and the maximum 
cardinality of 3 is thus 2 nK . We need a selection rule to select 2 nK ignition states among 2 n 
potential fine grains. Let 

K n -1 

3 = E ^2" (61) 

be the binary expansion of an integer j, (0 < j < 2 K — 1). For any fj, £ [0, K n — 1], there is a 
unique base- A" expansion of fi as 

n-l 

H = k i K * (0<h<K-l) 

i=0 

Define a selection rule a, 

a : [0, Kn - 11 [0, K" - 11 
(fe, i) — y fj, = o{k, i) 

A priori, any function a can be used. . . provided that the linear system, Eq. (|63p below, leads to 
feasible solutions. The simplest selection function is ai(k,i) = k + iK = 7. We have found by 
simulation that this simple function is sufficient, and more sophisticated selection rules only give 
a permutation of the label of the ignition states. Therefore, we will use this selection function 
throughout the rest of the paper, unless stated otherwise explicitly. 

Let j M be the coefficient of 2 M in the binary expansion of j in Eq. (|61[) . A convenient definition 
of the projection function is 

j) ( = f> 3v{k,i) = Jk+iK = jf (62) 
For ease of exposition, when no confusion is possible, we will write alternatively 

Now, we are going to show that this problem is roughly similar to a bipartite Bell-type system 
with Kn settings. 

Gauge probability distributions In the n-region-Tf-setting context, we will need nK gauge 
probability distributions g 7 (Aj) (or simply, g 7 (j)) for 7 = k + iK 6 [0, nK — 1] with j € [0,71 — 1] 
and k £ [0, K — 1]. Suppose that the observer in region lZi selects the setting fco- Let us derive 
the gauge distribution g 70 (j) for 70 = fco + ioK. Let p 70J be the unknown probability of Xj for 
j £ 3. Thus, we have to solve a linear system. 

For each setting vector u with u; = 6k and each outcome vector x, we will write one equation: 

Ui = 0k : Vui(i / io),Vx, : ^ Pioj = p ( x l u ) (63) 

jej(x,u) 

j G J(x, u) means that Vl : Xi — Hi(ui,\j), i.e., Xj £ (x|u). For each gauge distribution g 10 (j) 
we have card(3) unknowns p 10 j and 2 n K n ~ 1 equations among which it can be shown that only 
2(K + l) n_1 are independent. This system is generally degenerate and provides a set of solutions, 
S70C?) = Pioj- The computation fails if it is impossible to find non-negative solutions. Then, we 
have to try an indirect collapse, with more than one step. 

Consistency conditions The gauge distributions are not independent. As in Bell- type system, 
Eq. (|2U . we have a number of consistency conditions, translating complete local consistency, 
Eq, (l32|l . into gauge formulation. Whatever the gauge distribution g 7 with 7 = k + Ki, provided 
that Ui — k, we have: 

Prob(a;o; . . . ;av_i|u ; • • • ;Ur-i) = E &~<^ ( 64 ) 

\j&(x ;...:x r _ 1 \u \...;u r _ 1 ) 
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Measurement We proceed similarly as for totally correlated Bell-type pairs: In each region i, 
each observer selects his own setting m = 9k- Each configuration 74 = k + iK is transmitted with 
finite velocity towards the ignition point 3. Only one configuration, say 7i gn i, is kept at random. 
Now, at 3, we perform a trial in the ignition set A using the only probability distribution g 7igni (Aj) 
to draw a single ignition state Xj ; ni . Let j li be the coefficients of 2 7i in the binary expansion 
of Jigni- These coefficients, j li — il(7i,ji gn i) are transmitted with finite velocity towards the end 
regions IZi, where the final outcomes are thus Xi = II(7i, ji gn i). 

Conversely, let g7 igni (Aj) be any gauge distribution, with 7i gn j = fci gn i + ii sn iK. This gauge 
selection implies that the setting in region ii gn i is #fc igni . It is possible to compute the probability 
P(x|u). We have to collect all ignition states Aj within the target (x, u), given that ni = ni : 
Then, we have, 

P(x|u)= S7 lg »i(Ai) (65) 

where Aj £ (x|u) means that Wi,Xi = IL;(Mj, Aj). 



3.3.4 Multi-step collapse 

If we fail in computing a feasible set of gauge distributions, g 7 (Aj), we try a two-step collapse. 
This means a gauge selection of one region and a local trial to check again a new one-step collapse, 
but now in only n — I regions. If we fail again, we iterate for a three-step collapse, etc. Finally, 
any n-region system admits at least a classical collapse in n steps. 



4 Typical examples of multipartite systems 

We are now going to give some examplej"!. In bipartite systems, we will consider the cases of 
partially entangled regions as well as 'super-quantum' states. In tripartite system, we will deal 
with the GHZ and the W states. For simplicity, we will only compute these tripartite systems with 
two settings, namely, X and Y . We have already discussed of their entanglement behaviour with 
respect to the axis Z in Sec. 13.2.31 (Table [Til) . In each example, we will derive the multipartite 
entanglement entropies. 



4.1 One-region systems 

It may be convenient to have a general gauge formalism, valid in one region systems. Only the 



direct collapse is relevant. A one-region system with K settings (Table 14(a) I has K configurations 
and K gauges distributions. It requires a working set of 2 ignition states (Table [TT|) . as described 
in Table [T4(b)| 

(a) Probability distribution P(x|u). 
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(b) Gauge distributions gk(j) with 2 ignition states, j = and 
j = 2 K - 1 = 31. 
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Table 15: Typical one region system with K=5 settings. 



Such a system is trivially not entangled, and therefore the Bell-CHSH inequalities are irrelevant. 
Similarly, it is impossible to characterize a super-quantum behaviour, but a given probability 
distribution P(a;|'u) may force a number of constraints on the settings. Actually, when the system 
mimics a quantum situation, the targets (0\6k) and (l\9k) describe two qubits |£|#fc) with £ 6 {0, 1}, 
or alternatively two opposite unit vector Te^ on the surface of the Bloch sphere in R 3 , while the 
density operator of the system, p, is described by a vector r in the volume of the Bloch sphere still 
in R 3 . For instance, if po = pi — 1 in Table [T5l then the system is in the pure state \Q\0o) and 

In the previous versions of this paper, theses examples were widely erroneous. 
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|O|0i) and the constraint is that two settings do and 9\ have to be identical, i.e., eo = ei. More 
generally, the target probability is, 

P(O|0 fc ) = (1/2)(1 + r ■ e fc ) ; P(l]6» fc ) = (1/2)(1 - r • e fc ) 

where the dot stands for the ordinary scalar vector in R . Let r k be the projection of r on 

r k = r ■ e fc = 2P(O|0 fc ) -1 = 1- 2P(1|0&) 

Now, with a maximum of three independent settings, 6 kl , 8 k2 and 9k 3 , if we have 

ffej + rl 2 + rl 3 > 1 

the constraint is that the three vectors , efe 2 and , cannot be orthonormal in R 3 because the 
norme of r is less than 1. 

4.2 Bipartite systems 

We have already described Bell-type system within a particular formalism, devoted to totally 
correlated entities. Now, we will use the general formalism, valid for any degree of entanglement, 
including separable states. 

4.2.1 Separable states 

Consider the case of two independent entities £o and d. By definition, the entities are not 
entangled and the joint probability of the pair is defined by the product 

P(xo; xi|uo; «i) = Prob(a;o|ito) x Prob(xi|ui) 

The entangled entropy matrix of separable states is of course zero. 

The collapse mechanism involves physically two steps. Nevertheless, we have always found 
that a one-step collapse is successful. This opens the possibility of 'purification': a third region 
entangled with each of the two separable regions can be added to the system and a one-step overall 
collapse is likely to be triggered from this third region. 

4.2.2 General bipartite 2-setting systems 

Let K — 2. When we drop the condition of total correlation, Table [2] is expanded into Table \Wh , 
depending now on 8 parameters, gi, <J2, - • - , 98- It is possible to define a one-step collapse as follows: 
Each gauge distribution is defined by a linear system. We have 2x2 = 4 gauge distributions. For 
each gauge distribution, there are 16 unknowns for 6 independent equations when accounting for 
local consistency. We have then 16 ignition states but the rank of the linear system is 6 and thus 
each distribution can be computed with only 6 non zero components. It is easy to compute the 
partial distribution in each region (Tables [lllb and 116b ). The maximum bipartite entanglement 
entropy is equal to 1 bit. Especially, this maximum is obtained in some totally correlated systems 
like Bell's states but also in more exotic objects like the PR-Box (Table. 118b ) below. 



(a) P(xq;xi\uo;ui) 
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(b) Prob(x |«o)- 
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Table 16: General locally consistent 2-region system with 
to qg. All entries have to be non negative. 



(c) Prob(xi |ui). 
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2 possible settings Bq and 8i depending on 8 parameters qi 
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(a) Probabilities P(x|u) 
three settings (i 7^ j). 
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(b) Gauge distri- 
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(c) Entropy Matrix. 
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Table 17: Singlet state: Classical simulation of the two region singlet state with three settings 9q = X, d\ = Y and 
02 = Z: The six gauge distributions g-y(j) = g(j) are here identical. This means that the singlet state can be described 
by a quadruplet of equally likely 'hidden variables', namely A7, A28, A42, and A49. At each run, each observer submits 
its configuration labelled fi to a referee located at the ignition point. When she received the first configuration, the 
referee draws at random one of the four ignition states, say A Jigni , and transmits backwards to the observers the final 
outcomes xi = n(j'i gn i,7i). Such an isotropic pair of opposite spins in inconceivable in non-contextual systems. Finally, 
the singlet state is maximally entangled with a total entanglement of 1 bit. 



4.2.3 Revisiting Bell-type systems and the singlet state 

Of course, the general case includes the particular case of totally correlated Bell-type systems 
described in Sec. [2l and provides new gauge distributions in this context. This inflation of dis- 
tributions raises no real problem because ignition states and gauge distributions are not physical 
entities but only gauge parameters. For instance, we can recover the previous distributions of the 
EPR pairs (with the convention recalled in footnote [T]) by using a particular working set of 2K 
ignition states, namely j = (2 K + l)Di((r) (with < r < 2K — 1, Dk{t) being the double-plateau 
function of order K defined in Sec. I2.2.1[) . Then the gauge distributions gfe are duplicate and 

gk+K = gfc- 

Singlet state Specially, consider the singlet state, 

|VO = (1A/2)(|01> - |10», 

describing a pair of isotropic opposite spins with the usual convention, when e.g., |0) stands for 
spin up and 1 1} for spin down, irrespective of the region. This situation is generally considered as 
purely quantum and even inconceivable in classical physics. However, it is easily shown that this 
belief is baseless in contextual systems. Now, the probability distribution of the singlet state is 
given in Table [T7b for three settings, X, Y and Z. It is easy to compute a working set of gauge 
distributions. It is even possible to find a quadruplet of four ignition states, namely A7, A28, A42, 
and A49 using a single gauge distribution g(j) as shown in Table \T7b. Therefore, the singlet state 
is compatible with a set of four 'hidden variables' located at the ignition point. This proves that a 
potential pair of both isotropic and opposite vectors is surprisingly a classical although contextual 
concept. The total entanglement, Eq. (|50|) is maximum and equal to 1 bit, while the entropy matrix 
S2 is given in Table [T7b . Nevertheless, with the three settings X, Y and Z, the CHSH average is 
just equal to 2, and this entanglement is not detected by this criteriorP^I. 

4.2.4 PR-Box 

Another example of interest is the so-called 'PR-Box' proposed by Popescu and Rohrlich [3?3 | and 
recalled in Table H5b. This example is supposed to be hypothetical or 'super-quantum' because it 
exceeds the Tsirelson bound. Now, we see that the box corresponds to an 'ordinary' 2-region-2- 
setting-system, with qi — 53 = = q$ — q$ = 1/2, q2 ~ 1 and §5 = q*i = in Table [TU1 Define 
A = O , A' — 6\, B — 0i and B' = O in Eq.JS3] We have 

\S(A, B) + S(A', B) + S(A, B') - S{A' , B')\ = 4 > 2V2. 

Therefore, the system is 'super-quantum'. The collapse can be completed in classical physics, using 
e.g. the gauge distributions given in Table 118b . A possible implementation as a classical game 
using a fair coin is proposed in Table [181 caption. The system is clearly maximally entangled and 
its entanglement entropy is given in Table [T8b . This is an example of maximally entangled state 
which is not totally correlated (as defined in Sec. 12. 1 .3[) . The box is known to solve the problem 
of 'communication complexity' [67], [fj3 | in the sense that all distributed computations can be per- 
formed with a trivial amount of communication, i.e., with one bit [69]. This seems very surprising 
because the device can be easily implemented, but on the other hand, this classical realization 
involves a stage of classical communication. Finally, we have found that the box also describes 

However with four settings, it is well known that the Tsirelson bound 2\/2 is reached. 
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the second-step collapse of the 'super-quantum'-GHZ-states (Sec. 14.31 below), demonstrating the 
'super-quantum' behaviour of this state (see Sec. 13.3. 2[) . 



(a) P(x|u 



(b) Gauge distributions g-yfj) 
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(c) Entanglement entropy (in bits) 
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Table 18: PR-Box: (a) Probabilities P(x|u) proposed by Popescu and Rohrlich |3(I . This system is maximally 
entangled (c) and 'super-quantum'. The process can be implemented in classical physics, using the gauge distributions 
given in (b) with 4 ignition states (out of 16). This proves that the PR-box, while exceeding the Tsirelson bound 
Eo. 1 1531 . can be implemented in classical physics. We note that each distribution has only two equally likely outcomes. 
Therefore, it is possible to emulate the system as a simple game between two players Co and 0\ using a fair coin: Each 
player Oi chooses her/his setting ki = or 1 and sends to a referee (located at the ignition point) her/his configuration, 
7i = ki + 2i. The referee selects at random between 70 and 71 a gauge configuration, 7i gn i £ [0, 3], and tosses the coin. 
The coin is regarded as a gauge distribution g7 igni 0) with an assignment of the relevant ignition states j to 'head' and 
'tail' respectively, e.g., 'head' for j £ {3,6} and 'tail' for j £ {9, 15}. Then, the referee draws an ignition state j. For 
example, j = or 15 if 7i gn i = 0. At last, the final outcomes are Xi = Ii(~fi, j) respectively for i = 0,1. For example, if 
70 = 0, 71 = 1, 7i gn i = and j = 0, we have xq = x\ = 0. The entanglement entropy 82(^0! n i) = £( n 0i u l) ls given in 
(c). We have S2 = 1 = n — 1 where n = 2 is the number of regions. This proves that the system is maximally entangled 
while not totally correlated. 



4.3 Tripartite entangled systems 

It is well known that there are two irreducible families of tripartite quantum entangled systems, 
referred to as GHZ states and W-states. We will first analyse the W-states with only two settings. 

4.3.1 W-type states 

Consider an ensemble of three entities {£0, £1, £2}, placed respectively in three distant space 
regions 1Zq,1Z\ and IZ2 located in a plane {X, Y). The physical system is equivalent to a set of 
three spin 1/2 particles, initially entangled in the state 

= -U|001)z + |010) z + |100)z) 
V3 

specified along an axis Z perpendicular to the plane (X, Y), where, e.g., \Q)z stands for spin up 
and \l)z for spin down. Later the spins of the three particles are measured independently in the 
distant regions along either the axis X or Y. In each region, an observer selects freely a setting 
Ui {i = 0, 1 or 2). Therefore, we consider the tripartite system with only two settingj^l. Let 
= {X, Y}. The conditional probabilities P(a;o; x\; X2\uo\ ui; 112) = P(xju) are given in Table [T9b 
and clearly depends upon the settings in the three regions. Complete consistency holds and, e.g., 
the partial probability Prob(xo|i(o) in region TZq does not depend on the setting in 1Z\ and T&2- 
Thus, we have (Vwo, Vwi, Vtt2 G ©): 

1 1 1 
Prob(a;ojito) = ^ P(ko; xi;x2\uo; uv, 112) = - (66) 

and similarly for ports IZi and 1Z2- If we ignore the last port, the probability Prob(a;o; xi|uo; Ml) in 
regions (7?.o,7?.i), does not depend on the setting in 1Z2- Therefore, we have (\Mo , Vwi , Vti2 £ ©): 

1 

PTob{x ;xi\u Q ;u 1 ) = ^ P{x ;xi;X2\uo;ui;u 2 ) (67) 

The result is given in Table [20k . 
We have nevertheless considered the third setting Z in Sec. 13. 2.31 to compute its entanglement entropy (Table [T41 . 
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(a) W-state with two settings: Target probabilities P(x|u) 
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(b) Gauge distributions g 7 (j) 
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(c) Entanglement entropy (in bits) 
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Table 19: W-states: (a) Conditional probability P(xo; xv, X2\uq; ui; U2)- (b) Gauge distribution working set with 28 
ignition states Xj (out of 64) and 6 gauge distributions g 7 (j). (c) Tripartite entanglement entropy 83(1x0, ui , M2) and 
total entanglement £((uo, 1*1, 112)) (in bits). Since £(X, X, X) = £(Y, Y, Y) = 0.792 < 2 bits, this state is not maximally 
entangled. 
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(a) Partial probability Prob(xo; xi\uq; (b) Gauge probability P'(xq; xi\uq; «i). 
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(c) Partial entropy S2 (in bits) (d) Gauge entropy S' 2 (in bits) 
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Table 20: Partial and gauge bipartite subsystems in regions IZq and Tii of the W-state with two settings, (a) The partial 
probability in regions lio and TZi describes an entangled 2-region 2-setting system, (b) Gauge probability in regions IZq 
and Hi following a first collapse in region H2 with the setting 112 = X and the outcome X2 = 0. (c) The partial entropy 
computed from (a) shows that the maximum total entanglement is S2{X,X) = £(X, X) = S2{Y,Y) = £(Y, Y) = 0.35 
bit < 1 bit. The system is not maximally entangled, (d) The gauge entropy computed from (b) shows that the total 
entanglement is £' (X, X) = 0.09 bit while £'(Y, Y) = 0.35 bit is unchanged. 



This is the probability distribution of an entangled system obtained from the general 2-region 
2-setting system (Table fTok ) for gi = 52 = 7/12, 53 = q± = 3/4 and gs = ge = 96 = 97 = 98 = 1/12. 
It is easy to compute the bipartite entanglement entropy matrix. The result for one pair of regions 
is given in Table l2"Ub . As a result the bipartite degree of entanglement is e2 = 3 and the bipartite 
entanglement entropy S2, equal to the total entanglement is £(X,X) = £(Y,Y) = 0.35 < 1 bit. 
Therefore, the entanglement is not maximum and, for instance, is not detected by the CHSH 
inequalities. 

Two-step collapse It is possible to proceed to a partial collapse, e. g. in region 1Z2- We have 
to select a setting, e.g., u<z = X. The partial probability distribution in IZ2 is given by Eq. (|66[) 
and the two outcomes are equally likely. Suppose that we draw, e.g., the outcome 22 = 0. In 
regions IZq and 1Z\ the new probability distribution is easily computed and the result is given 
in Table [20b . This is a two-region two-setting system with gi = 1/4, g2 = qz — qi — 7/12 and 
9s = 96 = 97 = 98 = 1/12. We have computed also the total correlation in Table |2"UB . The collapse 
can be completed in one or two steps. 

One-step collapse Coming back to the full system, the number of configurations is nK = 
3x2 = 6. We can use the standard projection function, Eq. (|59[) . with just 2 nK — 2 6 = 64 
ignition states and 6 stochastic gauge distributions, g-y(Xj) with 7 = 0, . . . , 5. The result is given 
in Table [19b. 

4.3.2 GHZ-type states 

The second family of tripartite quantum entangled system is the GHZ-states, defined by Green- 
berger et al pFol ] . The physical system is a set of three spin 1/2 particles, initially entangled in the 
state 

|V) = 4= (1 000} z + |lll)z) (68) 
v2 

specifies along an axis Z perpendicular to the plane (X, Y). Again, the spins of the three particles 
are measured along either the axis X = 8q 01 the axis Y = 0\ and we only consider a 2-setting 
system, with © = {X, Y}. 

The target probabilities P(x|u) are given in Table [21k. In non-contextual systems, one would 
expect that each local outcome only depends on a local trial. This assumption is easily checked 
by assuming the existence of local random functions, defined as 3(lZi\9k) = x% when the entity £i 
is measured with the setting Oh- Let f(uo,ui,U2) = H(72,o|wo) + 3(7?-i|wi) + E(7?.2|it2) (mod 2). 
Clearly, f(X, X, X) = f(Y, Y, X) + f(Y, X, Y) + f(X, Y, Y) (mod 2). The assumption is falsified 
by Table [21k, since f(X,X,X) = 1 with certainty, while f(Y, Y, X) = 0, f(Y, X, Y) = and 
f(X, Y,Y) = with certainty. Therefore, the system is contextual. 

Complete local consistency holds. For instance, the partial probability Prob(a;o|uo) in region 
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(a) GHZ state with two settings: Target probabilities P(x|u) 
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(b) Gauge distributions g-y(j) 
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(c) Entropy (in bits) 
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Table 21: GHZ states for two settings X and Y : (a) Conditional probability P(xo; x\; X2\uq; ui; 112). Local realism 
is supposed to be falsified by computing a convenient function, /(«o, «i, M2) = a'o + xi + £'2 (mod 2). Whatever 
the outcomes, we see by inspection that f(X, X, X) = with certainty, while f(X,Y,Y) = 1, f(Y,X,Y) = 1 and 
f(Y,Y,X) = 1 with certainty, (b) Gauge distribution working set with 29 ignition states (out of 64). (c) Tripartite 
entanglement entropy Ss(uo, 111,112) and total entanglement £(«o, ^li "2) (in bits). Partial subsystems are separable. 
Therefore, bipartite entanglement entropies S2(ui,Uj) are identically zero. 



IZo does not depend on the setting in IZi and TI2 (Table [22k ). Thus, we have (Vlto, Vlti, V«2 £ ©): 

11 1 
Prob^ol^o) = ^2 ^2P(x ;xi;x2\u ;ui;u2) = - (69) 

and similarly for ports IZi and 1Z2- If we ignore the last port, the probability Prob(a;o; xi\uo;ui) in 
regions (72-0, 7£i), does not depend on the setting in 72.2- Therefore, we have (Vuo , Vlti , Vti2 G ©): 

1 1 
Prob(x ;xiju ;jii) = y] P(xo; xy, X2\uq; uy, tt a ) = 7 (70) 

This is obtained from the general 2-region 2-setting system (Table[l6k) for qi — q2 = qz — qi — 3/4 
and §6 = 96 = 97 = 98 = 1/4. Thus, the two regions are separable and the probabilities given by 
Eq (|70[) are simply the product of the probability obtained in Eq. (|69|l . Therefore the entanglement 
scheme is (0, 1). The bipartite entanglement entropy for any pair of region is zero (Table l22b ). By 
contrast, depending upon the settings, the signed tripartite entanglement entropy is equal to or 
— 1 bit and the total entanglement to or 1 bit(Table !21b ). We have seen that the GHZ-state is 
maximally entangled with the third setting (Z, Z, Z). 
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(a) Partial probability Prob(xo; xi\uo; 



(b) Gauge probability P' (xq; xi\uq; ui). 
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(c) Partial entropy S2 
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Table 22: Partial and gauge bipartite subsystems in regions TZq and "R.\ of the GHZ-state with two settings, (a) 
The partial probability in regions 7?.q and 7£i describes a separable 2-region 2-setting system, (b) Gauge probability in 
regions 1Zo and 7^i following a first collapse in region 7^2 with the setting U2 = X and the outcome X2 = 0. This is a 
Bell's state 0+ . (c) The partial subsystem (a) is separable. Therefore the partial entropy computed from (a) is zero, (d) 
By contrast, the gauge subsystem (b) is entangled. The gauge entropy (in bits) computed from (b) exhibits a maximum 
entanglement: S'{X,X) =£'(Y,Y) = 1 bit. 



Two-step collapse We can proceed to a 2-step collapse. Whatever the gauge region and the 
gauge setting, the bipartite gauge system is a Bell-state. Table [22b, d is computed for a collapse 
in region IZ2 with the setting 112 — X and the outcome X2 = 0. 

One-step collapse The one-step collapse rises no problem. The result is given in Table [2lb . 
4.3.3 Super-quantum GHZ state 

We are now going to describe a surprising 'super-quantum' GHZ-statJ^I. Actually, in the GHZ- 
state (Table [2lj) . the settings X and Y are not on an equal footing, even if this is not manifest in 
the definition, Eq. (|68[) . Now, if we demand a symmetrical behaviour between the two settings, 
we obtain the probability distribution of Table |2"3"R. that we will name Super- quantum GHZ state. 
It is easily shown that this distribution describes a locally consistent system. In addition, the 
partial probabilities are identical to the partial probability of the conventional GHZ-states and 
given by Eq. (|69l 17011 . As a result, the entanglement scheme is also (0,1). The interest of this 
example is threefold. Firstly, we will show that this state is indeed super-quantum by proceeding 
to a 2-step collapse. Secondly, this example cannot collapses in one step and this can be easily 
proved. Thirdly, we may construct a continuous set of similar modes, allowing an adjustment of 
the total entanglement in the range [0, 1] bit and a swiching between super-quantum, quantum 
and non contextual behaviours. 

We first try a one-step collapse. The number of configurations is 3 x 2 = 6. In order to classically 
encode Table [23k . we construct an ignition set A = with 2 6 = 64 ignition states, labelled 

j from j = to 63, and 6 stochastic gauge distributions, g 7 (Aj) with 7 = 0, . . . , 5. However, 
none of the six linear systems admits non negative solution. To circumvent this problem, we 
may try a different ignition set. The maximum number of distinct fine grains is 2 = 256. 
However, it is easily shown that these 256 fine grains are all included into at least one target of 
zero probability. Therefore, whatever the gauge, the probability of all ignition states will always 
be zero. In conclusion, within the realm of local realism, the super-quantum GHZ state cannot 
collapse in one step. 

We have thus to try a two-step collapse. Let X be the ignition point at the boundary of the 
three regions. Each observer selects freely his own setting, m = ki and sends his choice towards X, 
with a finite velocity. Suppose that the first received configuration is 112 = &2 from region 7?. 2 • As a 
result, we break the equilibrium at X from region 1Z2- We draw the outcome X2 = £2 from the local 
probability distribution Prob^lfe), Eq. (|69[l . The two possible outcomes £2 G {0, 1} are equally 
likely. Suppose for instance that £2 = 0. Next, in regions 1Zo,1Zi we compute the second-step 
probability distribution P'(xo; £i|tto; u i) — Prob(:ro; ari |^2 ; fe; ito; Mi)- The resulting distribution is 
given in Table l2"3"b for ft2 — and £2 = 0. On the other hand, if the observer O2 had chosen &2 = 1 
instead of &2 = 0, the distribution V' (xq;xi\uo;u\), now given in Table [23h for £2 = 0, would 
have been different. These probability distributions are identical to the probabilities of the PR- 
box, Table [TSk, (apart from the labelling). Therefore, this 'super-quantum' GHZ state is actually 

This example stems from a flaw in the previous versions of this paper. 
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(a) 'Super-quantum' GHZ state with two settings 9q and 9\: Probabilities P(x|u) 
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(b) Tripartite entropy S3 and total entanglement £ (in bits) 
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(c) Second-step probability P' when u 2 = 0q (d) Second-step probability P' when m 2 = 0i 
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(e) Gauge distributions for (c) (f) Gauge distributions for (d) 
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Table 23: 'Super-quantum' GHZ states: (a) Conditional probability P(xo; xv> X2 |mo; mi ; m 2 ) describing a 'super- 
quantum' GHZ system with two different settings 60 an d 01- It can be seen that this system is locally consistent 
and behaves exactly like the conventional GHZ state with respect to 'local realism', (b) Tripartite entanglement entropy 
S3 and total entanglement £ of the state. A one-step collapse is impossible. A partial collapse in region 7£ 2 leads to a 
PR-box. (c) Second-step probability P'(xo; xi|mo; mi) in regions 7?-o>7£i, following a gauge selection of region 7£ 2 and 
a free choice m 2 = 9q. (d) Second-step probability P'(a;o; x\ \uq\ u\) when the free choice is now m 2 = 0i. (e) Gauge 
distribution working set with 4 ignition states (out of 16) for the second-step probability P'(a;o; x\ |mo; mi) given in (c). 
This is a relabelling of Table [T8b . (f) Same as (e) for the distribution (d). 



'super-quantum'. To complete the two-step collapse, we have now to construct an ignition set for 
the two regions TZq and 72-1, accounting for the relevant distribution P'(a;o; a;i[«o; wi). This is easily 
achieved because the new system is actually a PR-box. Each gauge distribution requires just two 
ignition states, but we need four ignition states for the whole, as shown in Table l23e .f, This proves 
that the 'super-quantum' GHZ-states can be classically implemented, but only a two-step collapse 
is feasible. 

A continuous set of similar modes We have checked whether a quasi-super state can 
collapse in one step. Indeed, in the super-quantum GHZ state, a number of targets have zero 
probability (Table [2TV ). It is possible to replace by a small probability e as shown in Table [24] 
In spite of this change, the partial probabilities, Eq. (|69|) and (170|> are strictly conserved and 
therefore, any pair of regions remains separable. As a result, the entanglement scheme is still 
(0, 1), except when the three regions become independent. Thus, the new system may actually 
be viewed as a quasi-super GHZ state, provided that e remains small. By computer simulation, 
we have found that a one-step collapse is feasible only when e is in the range ]1/16,3/16[. For 
e = 2/16, the three regions become independent. 
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Table 24: 'Quasi-super-quantum' GHZ states: Conditional probability P(xo; xi; a^l^o; u 1; M2) of a quasi-GHZ system. 
Target probabilities are derived from the super-quantum GHZ state by replacing by e in Table l23t t. A one-step collapse 
is feasible for 1/16 < e < 3/16. Partial probabilities Eo J69l and 1701 are unchanged. 
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Behaviour of the 'quasi-super-quantum' object versus e 

For e = 1/16 = 0.0625, we have checked whether the system passes the GHZ test of 'local 
realism' (in fact, non-contextuality). With the same notations as above ( Table [2l"j) . the probability 
that f(9i, 0i, #1) = 1 (mod 2) is now 0.75 (instead of 1), while the probability that f(6o, do, 81) = 0, 
f(8o,0i,9o) = and / (#1, #o, #o) = is 0.56 (instead of 1). We can conclude that the behaviour of 
a quasi-GHZ state is similar to a strict mode, but there is no dramatic all-versus-nothing outcome 
as in the strict state. 

Finally, we have checked by simulation whether the system passes the Tsirelson test. We have 
found that a local collapse in one region gives a 'super-quantum' bipartite system when e < 0.0366. 
When e = 0.0366, the CHSH average (Eq. I53[l is about 2\/2. The total tripartite entanglement 
is then £ — 0.4 bit, irrespective of the settings. We can conclude that the system behaviour is 
'super-quantum' at least for < e < 0.0366. Note that this example indicates that a partial 
collapse can give an issue about the super-quantum behaviour of a tripartite system while direct 
3-region witnesses are unable to conclude. 

More generally, the current model should open the possibility to classically synthesize any 
locally consistent system, whether quantum or not, whatever its degree of entanglement. One 
has to define a number of regions, n, a number of settings, K, n local measurement entropies 
and C^)K m coefficients of entanglement for m = 2 to n. This is similar to the analysis of the 
so-called 'nonsignaling polytope' in quantum information. This discussion is beyond the scope of 
the present paper. 

5 Conclusion 

Contextuality is not the privilege of the quantum world. Quite the reverse, we have shown in a 
previous paper [2(| that a strictly classical dice game can be contextualfy dependent and even 
nonsignaling in the quantum sense. Therefore, the difficult task to reconcile quantum mechanics 
with non-contextual logic is unnecessary. In the current paper we have constructed a model named 
'stochastic gauge system' to compute contextually dependent classical systems. We have shown 
that the theory described quantum randomness as well and can be used to simulate quantum 
states and even non-local boxes in classical physics, e.g., EPR pairs, GHZ states or PR-boxes. 
Quantum systems are compared with classical extended objects in equilibrium and quantum col- 
lapses are therefore identified with classical breaks of equilibrium. Furthermore, the model gives 
a straightforward classical interpretation of entanglement and entanglement entropy and provides 
simple tools for characterizing multipartite entanglements. 
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In our opinion, the present theory opens the way for a number of important questionings, 
regarding quantum theory, quantum information and even the general conception of physics, ob- 
viously beyond the scope of this paper. 

In quantum mechanics, we have shown that classical local realism is compatible with the 
conventional formulation of the theory, contrary to a general belief. This result is likely to explain 
a number of difficulties and opens up also new topics. For instance, in the relativistic domain, 
the conventional interpretation fails to define a consistent probability distribution. In our opinion, 
the concept of gauge probability should be a key ingredient to settle this problem. The notion 
of quantum collapse should be also revisited and compared with a classical break of equilibrium. 
New issues, like collapse kinetics should probably be investigated. Incidentally, while a verification 
of Bell's inequality violation is useless, a meaningful test to check the present theory would be 
an experimental probe of the collapse kinetics between entangled regions, in terms, e.g., of the 
distance between particles. 

The model proposes also an interpretation of quantum parallelism. A quantum system can be 
viewed as a complex entangled object in equilibrium. When the equilibrium is broken, a number 
of new situations become potentially possible. Quantum parallelism describes the set of all these 
potential situations. In quantum information theory, a final measurement will draw just one 
outcome. The theory uses entanglement as a resource, but again, this is not a quantum privilege. 
In our opinion, contextuality should be a resource in classical information as well. Indeed, 'quantum 
randomness' can be viewed as strictly classical, provided that one includes an additional concept in 
the conventional stochastic theory, namely, the notion of gauge distribution. Classical simulation 
of quantum algorithms should open up new paths in classical computation, beside deterministic 
and probabilistic algorithms. More generally, a number of notions like quantum cryptography, 
quantum complexity or communication complexity (see above Sec. 14. 2p should probably deserve 
to be revisited. 

Nowadays, the principle of gauge invariance plays a key role in modern physics. This is not 
very surprising because every physical measurement is performed within a particular framework 
while physical reality has to be independent of any referential and any system of units. Therefore, 
mathematical tools have been devised: Dimensional analysis, covariance principle, and more gen- 
erally, gauge theories. The gauge invariance principle has proved to be extraordinary powerful in 
all branches of physics. 

Surprisingly, the measurement in quantum theory remains nevertheless thought as absolute. 
The seed of this conception lies probably in the very origin of the theory, i.e., the Copenhagen 
interpretation of quantum mechanics, according to which the squared amplitude of the wave func- 
tion is an absolute probability. In our opinion, this belief is inconsistent, as highlighted by the EPR 
paradox or the problematic generalization of this interpretation in special relativity. 

The view of an absolute probability space is also the origin of a controversy between 'Bayesians' 
and 'Frequentists' in the probability domain. In our opinion, the 'Frequentist' position is simply not 
sustainable. For instance, we have emphasized in the present paper that a number of fundamental 
concepts like partial measurements are subjective by definition. 

Actually, physics aims to understand the entirety of the universe, in which a number of long 
range interactions are recorded. By contrast, each observer is embedded into his own limited 
causal bubble or causal diamond, bounded by light sheets, black hole and cosmological horizons. 
Therefore, he can only describe a finite part of the world and, apart from long range interactions, 
he has to account for the rest of the universe by his best evaluation of his irreducible ignorance, 
technically in form of entropy located on his causal horizons. 

This vision of the universe seems very new in physics: Unpredictability is no more in contra- 
diction with determinism but, quite the reverse, an inevitable consequence. This implies a break 
between the general Laplace concept of determinism and the technical notion of causality. In this 
respect, the dispute concerning the completeness of quantum mechanics appears totally unrealistic. 
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